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ABSTRACT 

The  infrared  structure  of  quantum  gravity  is  explored  by  solving  a  lattice  version  of  the  Wheeler- 
DeWitt  equations.  In  the  present  paper  only  the  case  of  2  + 1  dimensions  is  considered.  The  nature 
of  the  wave  function  solutions  is  such  that  a  finite  correlation  length  emerges  and  naturally  cuts 
off  any  infrared  divergences.  Properties  of  the  lattice  vacuum  are  consistent  with  the  existence 
of  an  ultraviolet  fixed  point  in  G  located  at  the  origin,  thus  precluding  the  existence  of  a  weak 
coupling  perturbative  phase.  The  correlation  length  exponent  is  determined  exactly  and  found 
to  be  v  =  6/11.  The  results  obtained  so  far  lend  support  to  the  claim  that  the  Lorentzian  and 
Euclidean  formulations  belong  to  the  same  field-theoretic  universality  class. 
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1  Introduction 


It  is  possible  that  the  well-known  ultraviolet  divergences  affecting  the  perturbative  treatment  of 
quantum  gravity  in  four  dimensions  point  to  a  fundamental  vacuum  instability  of  the  full  theory. 
If  this  is  the  case,  then  the  correct  identification  of  the  true  ground  state  for  gravitation  necessarily 
requires  the  introduction  of  a  consistent  nonperturbative  cutoff.  To  this  day  the  only  known  way 
to  do  this  reliably  in  quantum  field  theory  is  via  the  lattice  formulation.  Nevertheless,  previous 
work  on  lattice  quantum  gravity  has  dealt  almost  exclusively  with  the  Euclidean  formulation  in 
d  dimensions,  treated  via  the  manifestly  covariant  Feynman  path  integral  method.  Indeed  the 
latter  is  very  well  suited  for  numerical  integration,  and  many  analytical  and  numerical  results  have 
been  obtained  over  the  years  within  this  framework.  However  the  issue  of  their  relationship  to 
the  Lorentzian  theory  has  remained  largely  open,  at  least  from  the  point  of  view  of  a  rigorous 
treatment.  The  main  supporting  arguments  for  the  Euclidean  approach  come  from  the  fact  that 
the  above  equivalence  holds  true  for  other  field  theories  (no  exceptions  are  known),  and  from  the 
fact  that  in  gravity  itself  it  is  rigorously  true  to  all  orders  in  the  weak  field  expansion. 

In  this  paper  we  will  focus  on  the  Hamiltonian  approach  to  gravity,  which  assumes  from  the 
beginning  a  metric  with  Lorentzian  signature.  In  order  to  obtain  useful  insights  regarding  the  non¬ 
perturbative  ground  state,  a  Hamiltonian  lattice  formulation  was  introduced  based  on  the  Wheeler- 
DeWitt  equation,  where  the  quantum  gravity  Hamiltonian  is  written  down  in  the  position-space 
representation.  In  a  previous  paper  [1]  a  general  discrete  Wheeler-DeWitt  equation  was  given 
for  pure  gravity,  based  on  the  simplicial  lattice  formulation  originally  developed  by  Regge  and 
Wheeler.  On  the  lattice  the  infinite-dimensional  manifold  of  continuum  geometries  is  replaced  by 
a  finite  manifold  of  piecewise  linear  spaces,  with  solutions  to  the  lattice  equations  then  providing 
a  suitable  approximation  to  the  continuum  gravitational  wave  functional.  The  lattice  equations 
were  found  to  be  explicit  enough  to  allow  the  development  of  potentially  useful  practical  solutions. 
As  a  result,  a  number  of  sample  quantum  gravity  calculations  were  carried  out  in  2  +  1  and  3  +  1 
dimensions.  These  were  based  mainly  on  the  strong  coupling  expansion  and  on  the  Rayleigh- Ritz 
variational  method,  the  latter  implemented  using  a  set  of  correlated  product  (Slater-Jastrow)  wave 
functions. 

Here,  we  extend  the  work  initiated  in  [1]  and  show  how  exact  solutions  to  the  lattice  Wheeler- 
DeWitt  equations  can  be  obtained  in  2  +  1  dimensions  for  arbitrary  values  of  Newton’s  constant 
G.  The  procedure  we  follow  is  to  solve  the  lattice  equations  exactly  for  several  finite  regular 
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triangulations  of  the  sphere  and  then  extend  the  result  to  an  arbitrarily  large  number  of  triangles. 
One  finds  that  for  large  enough  areas  the  exact  lattice  wave  functional  depends  on  geometric 
quantities  only,  such  as  the  total  area  and  the  total  integrated  curvature  (which  in  2  +  1  dimensions 
is  just  proportional  to  the  Euler  characteristic).  The  regularity  condition  on  the  solutions  of  the 
wave  equation  at  small  areas  is  shown  to  play  an  essential  role  in  constraining  the  form  of  the  wave 
functional,  which  we  eventually  find  to  be  expressible  in  closed  form  as  a  confluent  hyper  geometric 
function  of  the  first  kind.  Later  it  will  be  shown  that  the  resulting  wave  function  allows  an  exact 
evaluation  of  a  number  of  useful  (and  manifestly  diffeomorphism-invariant)  averages,  such  as  the 
average  area  of  the  manifold  and  its  fluctuation. 

From  these  results  a  number  of  suggestive  physical  results  can  be  obtained,  the  first  one  of  which 
is  that  the  correlation  length  in  units  of  the  lattice  spacing  is  found  to  be  finite  for  all  G  >  0,  and 
diverges  at  G  =  0.  Such  a  result  can  be  viewed  as  consistent  with  the  existence  of  an  ultraviolet 
fixed  point  (or  a  phase  transition  in  statistical  field  theory  language)  in  G  located  at  the  origin,  thus 
entirely  precluding  the  existence  of  a  weak  coupling  phase  for  gravity  in  2  +  1  dimensions.  Simple 
renormalization  group  arguments  would  then  suggest  that  gravitational  screening  is  not  physically 
possible  in  2  +  1  dimensions,  and  that  gravitational  antiscreening  is  the  only  physically  realized 
option  in  this  model.  A  second  result  that  follows  from  our  analysis  is  an  exact  determination  of  the 
critical  correlation  length  exponent  for  gravity  in  2  +  1  dimensions,  which  is  found  to  be  v  =  6/11. 
It  is  known  that  the  latter  determines,  through  standard  renormalization  group  arguments,  the 
scale  dependence  of  the  gravitational  coupling  in  the  vicinity  of  the  ultraviolet  fixed  point. 

A  short  outline  of  the  paper  is  as  follows.  In  Sec.  2,  as  a  general  background  to  the  rest  of 
the  paper,  we  briefly  describe  the  formalism  of  classical  canonical  gravity,  as  originally  formulated 
by  Arnowitt,  Deser  and  Misner.  The  continuum  Wheeler-DeWitt  equation  and  its  invariance 
properties  are  introduced  as  well  at  this  stage.  In  Sec.  3  we  introduce  the  lattice  Wheeler-DeWitt 
equation  derived  in  a  previous  paper  |T],  and  later  Sec.  4  makes  more  explicit  various  quantities 
appearing  in  it.  This  last  section  also  discusses  briefly  the  role  of  continuous  lattice  diffeomorphism 
invariance  in  the  Regge  framework  as  it  applies  to  the  present  case  of  2  +  1-dimensional  gravity.  Sec. 
5  focuses  on  the  scaling  properties  of  the  lattice  equations  and  various  sensible  choices  for  the  lattice 
coupling  constants,  with  the  aim  of  giving  eventually  a  more  transparent  form  to  the  wave  function 
results.  Sec.  6  gives  a  detailed  outline  of  the  general  method  of  solution  for  the  lattice  equations  and 
then  gives  the  explicit  solution  for  a  number  of  regular  triangulations  of  the  sphere.  Later,  a  general 
form  of  the  wave  function  is  given  that  covers  all  the  previous  discrete  cases  and  allows  a  subsequent 
study  of  the  infinite  volume  limit.  Sec.  7  focuses  on  one  of  the  simplest  diffeomorphism-invariant 
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averages  that  can  be  computed  from  the  wave  function,  namely  the  average  total  area.  A  brief 
discussion  follows  on  how  the  latter  quantity  relates  to  the  corresponding  averages  computed  in  the 
Euclidean  theory.  Sec.  8  extends  the  calculation  to  the  area  fluctuation  and  shows  how  the  critical 
exponents  (anomalous  dimensions)  of  the  2  +  1-gravity  theory  can  be  obtained  from  the  exact  wave 
function  solution,  using  some  rather  straightforward  scaling  arguments.  Sec.  9  discusses  some 
simple  physical  implications  that  can  be  inferred  from  the  values  of  the  exact  exponents  and  the 
fact  that  quantum  gravity  in  2  +  1  dimensions  does  not  seemingly  possess,  in  either  the  Euclidean  or 
Lorentzian  formulation,  a  weak  coupling  phase.  Sec.  10  contains  a  summary  of  the  results  obtained 
so  far. 


2  Continuum  Wheeler-DeWitt  Equation 


Since  this  paper  involves  the  canonical  quantization  of  gravity  we  begin  here  with  a  very  brief 
summary  of  the  classical  canonical  formalism  [2]  as  derived  by  Arnowitt,  Deser  and  Misner  [3]. 
While  many  of  the  results  presented  in  this  section  are  rather  well  known,  it  will  be  useful,  in  view 
of  later  applications,  to  recall  the  main  results  and  formulas  and  provide  suitable  references  for 
expressions  used  later  in  the  paper. 

The  first  step  in  developing  a  canonical  formulation  for  gravity  is  to  introduce  a  time-slicing 
of  space-time,  by  introducing  a  sequence  of  spacelike  hypersurfaces  labeled  by  a  continuous  time 
coordinate  t.  The  invariant  distance  is  then  written  as 

ds2  =  — dr 2  =  g^dx^dx"  =  gij  dxl  dx^  +  2 gijNldx^dt  —  ( N 2  —  gij  NqN^)dt2  ,  (1) 


where  xl  ( i  =  1,2,3)  are  coordinates  on  a  three-dimensional  manifold  and  r  is  the  proper  time,  in 
units  with  c  =  1 . 

Indices  are  raised  and  lowered  with  gjj (x)  (• i,j  =  1,  2, 3),  which  denotes  the  three- metric  on  the 
given  spacelike  hypersurface,  and  N (x)  and  Nl(x)  are  the  lapse  and  shift  functions,  respectively. 
It  is  customary  to  mark  four-dimensional  quantities  by  the  prefix  4  so  that  all  unmarked  quantities 
will  refer  to  three  dimensions  (and  are  occasionally  marked  explicitly  by  a  3).  In  terms  of  the 
original  four-dimensional  metric  Ag^u  one  has 


% o  '//oA  _  ( NkNk  —  N2  Nj 
'.+o  l!Jij 



Ni  gtj  j 

which  then  gives  for  the  spatial  metric  and  the  lapse  and  shift  functions 


9ij  —  9ij 


4  jxn-1/2 


n  =  (-y°) 


Ni  =  4g0i  . 


(2) 


(3) 
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For  the  volume  element  one  has 


V~4g  =  N  Va 


(4) 


where  the  latter  involves  the  determinant  of  the  three- metric,  g  =  detgij.  As  usual  glJ  denotes  the 
inverse  of  the  matrix  gtJ . 

A  transition  from  the  classical  to  the  quantum  description  of  gravity  is  obtained  by  promoting 
the  metric  gij,  the  conjugate  momenta  it13 ,  the  Hamiltonian  density  H  and  the  momentum  density 
Hi  to  quantum  operators,  with  g ^  and  nlJ  satisfying  canonical  commutation  relations.  In  particular, 
the  classical  constraints  now  select  a  physical  vacuum  state  |  ,  such  that  in  the  source- free  case 


H\V)  =  0  Hi  I’F)  =  0 


(5) 


and  in  the  presence  of  sources  more  generally 

T|\H)  =  0  Ti  |*>  =  0  ,  (6) 

where  T  and  Tj  now  include  matter  contributions  that  should  be  added  to  H  and  Hi.  The  momen¬ 
tum  constraint  involving  Hi  or  more  generally  Ti,  ensures  that  the  state  functional  does  not  change 
under  a  transformation  of  coordinates  xl ,  so  that  T  depends  only  on  the  intrinsic  geometry  of  the 
3-space.  The  Hamiltonian  constraint  is  then  the  only  remaining  condition  that  the  state  functional 
must  satisfy. 

As  in  ordinary  nonrelativistic  quantum  mechanics,  one  can  choose  different  representations  for 
the  canonically  conjugate  operators  g^  and  n13 .  In  the  functional  position  representation  one  sets 

9ij(x)  ->•  ^(x)  7 ru(x)  ->■  -ih  •  16? tG  •  §  .  (7) 

OQij  lxj 

In  this  picture  quantum  states  become  wave  functionals  of  the  three- metric  g,j(x), 

l^>  ^  bij(x)]  ■  (8) 

The  two  quantum-constraint  equations  in  Eq.  ([6|)  then  become  the  Wheeler-DeWitt  equation  [3J  E 

m 

{-16lra'Gii'k,s£^l  ~  -  2A>  +  **}  4'fe(x)1  =  °'  (9) 

and  the  momentum  constraint  listed  below.  Here  is  the  inverse  of  the  DeWitt  supermetric, 

given  by 

Gij^ki  —  2  9  /  ( 9ik9ji  T  9u9jk  T  a>gijgki)  ,  (lb) 
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with  parameter  a  =  —  1.  The  three-dimensional  version  of  the  DeWitt  supermetric  itself,  G^,kl{x) 
is  given  by 

GljM  =  \Va  (glkgjl  +  gagJk  +  ugljgkl)  ,  (n) 

with  parameter  a  in  Eq.  (1101)  related  to  a  in  Eq.  m  by  a  =  — 2a/(2  +  3a),  so  that  a  =  — 1  gives 
d  =  —2  (note  that  this  is  dimension  dependent).  In  the  position  representation  the  diffeomorphism 
(or  momentum)  constraint  reads 

1 2  i  9ij  Vk  +  Hf  |  *[gij (x)]  =  0  ,  (12) 

where  H ^  and  Hf  are  possible  matter  contributions.  In  the  following,  we  shall  set  both  of  these  to 
zero,  as  we  will  focus  here  almost  exclusively  on  the  pure  gravitational  case. 

A  number  of  basic  issues  need  to  be  addressed  before  one  can  gain  a  full  and  consistent  un¬ 
derstanding  of  the  dynamical  content  of  the  theory  (see,  for  example,  m  0  n  na  mi  as  a  small 
set  of  representative  references).  These  include  possible  problems  of  operator  ordering,  and  the 
specification  of  a  suitable  Hilbert  space,  which  entails  at  some  point  a  choice  for  the  inner  product 
of  wave  functionals,  for  example  in  the  Schrodinger  form 

(T|<f>)  =  J  d/4g]  <I>fe]  (13) 

where  dfi[g\  is  some  appropriate  measure  over  the  three-metric  g.  Note  also  that  the  continuum 
Wheeler-DeWitt  equation  contains,  in  the  kinetic  term,  products  of  functional  differential  oper¬ 
ators  which  are  evaluated  at  the  same  spatial  point  x.  One  would  expect  that  such  terms  could 
produce  6 ®  (0)  -type  singularities  when  acting  on  the  wave  functional,  which  would  then  have  to  be 
regularized  in  some  way.  The  lattice  cutoff  discussed  below  is  one  way  to  provide  such  an  explicit 
ultraviolet  regularization. 

A  peculiar  property  of  the  Wheeler-DeWitt  equation,  which  distinguishes  it  from  the  usual 
Schrodinger  equation  H^>  =  ihdt'f’,  is  the  absence  of  an  explicit  time  coordinate.  As  a  result,  the 
rhs  term  of  the  Schrodinger  equation  is  here  entirely  absent.  The  reason  is  of  course  diffeomorphism 
invariance  of  the  underlying  theory,  which  expresses  now  the  fundamental  quantum  equations  in 
terms  of  fields  and  not  coordinates. 


3  Lattice  Hamiltonian  for  Quantum  Gravity 

In  constructing  a  discrete  Hamiltonian  for  gravity,  one  has  to  decide  first  what  degrees  of  freedom 
one  should  retain  on  the  lattice.  One  possibility,  which  is  the  one  we  choose  to  pursue  here,  is  to 
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use  the  more  economical  (and  geometric)  Regge- Wheeler  lattice  discretization  for  gravity  [12,  13] . 
with  edge  lengths  suitably  defined  on  a  random  lattice  as  the  primary  dynamical  variables.  Even 
in  this  specific  case  several  avenues  for  discretization  are  possible.  One  could  discretize  the  theory 
from  the  very  beginning,  while  it  is  still  formulated  in  terms  of  an  action,  and  introduce  for  it  a 
lapse  and  a  shift  function,  extrinsic  and  intrinsic  discrete  curvatures  etc.  Alternatively  one  could 
try  to  discretize  the  continuum  Wheeler-DeWitt  equation  directly,  a  procedure  that  makes  sense  in 
the  lattice  formulation,  as  these  equations  are  still  given  in  terms  of  geometric  objects,  for  which 
the  Regge  theory  is  very  well  suited.  It  is  the  latter  approach  which  we  will  proceed  to  outline  here. 

The  starting  point  for  the  following  discussion  is  therefore  the  Wheeler-DeWitt  equation  for 
pure  gravity  in  the  absence  of  matter,  Eq.  ([9]), 

{  “ (16,rG,)2 G«-ii(x)  ^  <  3r{x)  “  2A >  }  'f[9«(x)1  =  0  <14> 
and  the  diffeomorphism  constraint  of  Eq.  m, 

|  2  tgij(x)  Vfc(x)  ^  ^x)  |  ^bu(x)]  =  0.  (15) 

Note  that  these  equations  express  a  constraint  on  the  state  |T)  at  every  x,  each  of  the  form 
tf(x)  |T)  =  0  and  H%  (x)|T)  =  0. 

On  a  simplicial  lattice  M  ES  US  US  HE]  (see  for  example  m,  and  references  therein,  for  a 
more  complete  discussion  of  the  lattice  formulation  for  gravity)  one  knows  that  deformations  of  the 
squared  edge  lengths  are  linearly  related  to  deformations  of  the  induced  metric.  In  a  given  simplex 
(7,  take  coordinates  based  at  a  vertex  0,  with  axes  along  the  edges  from  0.  The  other  vertices  are 
each  at  unit  coordinate  distance  from  0  (see  Figs.  1,  2  and  3  as  an  example  of  this  labeling  for  a 
triangle).  In  terms  of  these  coordinates,  the  metric  within  the  simplex  is  given  by 

l  {loi  +  lOj  ~lij)  ■  (16) 

Note  that  in  the  following  discussion  only  edges  and  volumes  along  the  spatial  direction  are  involved. 
It  follows  that  one  can  introduce  in  a  natural  way  a  lattice  analog  of  the  DeWitt  supermetric  of 
Eq.  ([Ill)  by  adhering  to  the  following  procedure  [201  [21j.  First  one  writes  for  the  supermetric  in 
edge  length  space 

||  SI2  ||2  =  Gij{l2)  M2  M)  ,  (17) 

ij 

with  the  quantity  Gl^(l2)  suitably  defined  on  the  space  of  squared  edge  lengths.  By  a  straightforward 
exercise  of  varying  the  squared  volume  of  a  given  simplex  a  in  d  dimensions 

v2(° )  =  (ar)2det^V))  (18) 
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to  quadratic  order  in  the  metric  (on  the  rhs),  or  in  the  squared  edge  lengths  belonging  to  that 
simplex  (on  the  lhs),  one  is  led  to  the  identification 


Gij{l2) 


d\  y  — —  92  v 

■  4"  V(a)  dlf  81] 


(19) 


It  should  be  noted  that  in  spite  of  the  appearance  of  a  sum  over  simplices  a,  GlHj2)  is  local,  since 
the  sum  over  a  only  extends  over  those  simplices  which  contain  either  the  i  or  the  j  edge. 

At  this  point  one  is  finally  ready  to  write  a  lattice  analog  of  the  Wheeler-DeWitt  equation  for 
pure  gravity,  which  reads 

|  —  (167tG)2  Gij(l2)  ~  'M1)  t'W2)  -  2A]  |  ni2}  =  0  ,  (20) 

with  Gij{l2)  the  inverse  of  the  matrix  G2i (l2)  given  above.  The  range  of  the  summation  over  i  and 
j  and  the  appropriate  expression  for  the  scalar  curvature,  in  this  equation,  are  discussed  below  and 
made  explicit  in  Eq.  (ED- 

Equations  ([9])  or  (l20l)  express  a  constraint  equation  at  each  “point”  in  space.  Here  we  will 
elaborate  a  bit  more  on  this  point.  On  the  lattice,  points  in  space  are  replaced  by  a  set  of  edge 
labels  i,  with  a  few  edges  clustered  around  each  vertex  in  a  way  that  depends  on  the  dimensionality 
and  the  local  lattice  coordination  number.  To  be  more  specific,  the  first  term  in  Eq.  (1201)  contains 
derivatives  with  respect  to  edges  i  and  j  connected  by  a  matrix  element  Gjj  which  is  nonzero  only 
if  i  and  j  are  close  to  each  other,  essentially  nearest  neighbor.  One  would  therefore  expect  that  the 
first  term  could  be  represented  by  just  a  sum  of  edge  contributions,  all  from  within  one  (d  —  1)- 
simplex  <7  (a  tetrahedron  in  three  dimensions).  The  second  term  containing  3R(l2)  in  Eq.  (1201)  is  also 
local  in  the  edge  lengths:  it  only  involves  a  handful  of  edge  lengths,  which  enter  into  the  definition 
of  areas,  volumes  and  angles  around  the  point  x,  and  follows  from  the  fact  that  the  local  curvature 
at  the  original  point  x  is  completely  determined  by  the  values  of  the  edge  lengths  clustered  around 
i  and  j.  Apart  from  some  geometric  factors,  it  describes,  through  a  deficit  angle  5h,  the  parallel 
transport  of  a  vector  around  an  elementary  dual  lattice  loop.  It  should,  therefore,  be  adequate  to 
represent  this  second  term  by  a  sum  over  contributions  over  all  (d  —  3)-dimensional  hinges  (edges 
in  3+1  dimensions)  h  attached  to  the  simplex  a,  giving,  therefore,  in  three  dimensions 


(16vrG)2  Y,  GH  (*) 


d2 


i,jCa 


di2di2 


2  nah  lh  3h  +  2A  Va 

hCcr 


Tin  =  0 


(21) 


Here  5h  is  the  deficit  angle  at  the  hinge  h,  lh  the  corresponding  edge  length,  and  Va  =  yf  g(a)  the 
volume  of  the  simplex  (tetrahedron  in  three  spatial  dimensions)  labeled  by  a.  Gij  (a)  is  obtained 
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either  from  Eq.  (1191)  or  from  the  lattice  transcription  of  Eq.  Col) 

Gij,kl(a)  =  ^9~L/2(cr)[9ik{cr)gji(cr)  +  ga(a)gjk(a)  -  gij(a)gki(cr)\  ,  (22) 

with  the  induced  metric  gij(cr)  within  a  simplex  a  given  in  Eq.  (1161).  The  combinatorial  factor  nak 
ensures  the  correct  normalization  for  the  curvature  term,  since  the  latter  has  to  give  the  lattice 
version  of  /  sfg^R  =  2  Yh  (hi  three  spatial  dimensions)  when  summed  over  all  simplices  a.  The 
inverse  of  nak  counts,  therefore,  the  number  of  times  the  same  hinge  appears  in  various  neighboring 
simplices  and  consequently  depends  on  the  specific  choice  of  underlying  lattice  structure;  for  a  flat 
lattice  of  equilateral  triangles  in  two  dimensions,  nak  =  l/6@  The  lattice  Wheeler-DeWitt  equation 
given  in  Eq.  (1211)  was  the  main  result  of  a  previous  paper  |T] . 


4  Explicit  Setup  for  the  Lattice  Wheeler-DeWitt  Equation 

In  this  section,  we  shall  establish  our  notation  and  derive  the  relevant  terms  in  the  discrete  Wheeler- 
DeWitt  equation  for  a  simplex.  From  now  on  we  shall  focus  almost  exclusively  on  the  case  of  2+1 
dimensions.  The  basic  simplex  in  this  case  is,  of  course,  a  triangle,  with  vertices  and  squared  edge 
lengths  labelled  as  in  Fig.  1.  We  set  l qX  =  a,  l\2  =  b ,  Iq2  =  c-  The  components  of  the  metric  for 
coordinates  based  at  vertex  0,  with  axes  along  the  01  and  02  edges,  are 


5ii  =  a,  912  =  a  +  c-b ),  g22  =  c. 


(23) 


The  area  A  of  the  triangle  is  given  by 


A2  =  -}-[  2  (ab  +  bc  +  ca)  —  a2  —  b2  —  c2]  , 

16 


(24) 


so  the  supermetric  Gli ,  according  to  Eq.  (USD,  is 


1 


1  -1  -V 


Gl>  =  —  -i  1  -1 

4A  l-i  -i  i 


(25) 


Thus  for  the  triangle  we  have 


G 


d2 


lJ  dsi  dsj 


=  -4  A 


d2 


+ 


d 2 


da  db  db  dc 


+ 


d2 


dc  da 


and  the  Wheeler-DeWitt  equation  is 


(26) 


(16? tG)2  4 A 


d2 


+ 


d2 


+ 


d2 


da  db  db  dc  dc  da 


2  Tl (jh  J2Sh  +  2A4U[s]  =  0,  (27) 
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2 


Figure  1:  A  triangle  with  labels. 


Figure  2:  Neighbors  of  a  given  triangle.  The  picture  illustrates  the  fact  that  the  Laplacian  A (72) 
appearing  in  the  kinetic  term  of  the  lattice  Wheeler-DeWitt  equation  (here  in  2+1  dimensions) 
contains  edges  a,  b ,  c  that  belong  both  to  the  triangle  in  question,  as  well  as  to  several  neighboring 
triangles  (here  three  of  them)  with  squared  edges  denoted  sequentially  by  s\  =  l\  ...  s§  = 
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Figure  3:  A  small  section  of  a  suitable  dynamical  spatial  lattice  for  quantum  gravity  in  2+1 
dimensions. 
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where  the  sum  is  over  the  three  vertices  h  of  the  triangle. 

In  the  following  sections  we  will  be  concerned  at  some  point  with  various  discrete,  but  generally 
regular,  triangulations  of  the  two-sphere,  such  as  the  tetrahedron,  the  octahedron  and  the  icosahe¬ 
dron.  These  were  already  studied  in  some  detail  in  [221  [23 j.  A  key  aspect  of  the  Regge  theory  is 
the  presence  of  a  continuous,  local  lattice  diffeomorphism  invariance,  whose  main  aspects  in  regard 
to  their  relevance  for  the  3  +  1  formulation  of  gravity  were  already  addressed  in  some  detail  in 
[1J  in  the  context  of  the  lattice  weak  field  expansion.  Here  we  will  add  some  remarks  about  how 
this  local  invariance  manifests  itself  in  the  2  +  1  formulation,  and,  in  particular,  for  the  discrete 
triangulations  of  the  sphere  studied  later  on  in  this  paper.  Of  some  relevance  is  the  presence  of 
exact  zero  modes  of  the  gravitational  lattice  action,  reflecting  a  local  lattice  diffeomorphism  invari¬ 
ance,  present  already  on  a  finite  lattice.  Since  the  Einstein  action  is  a  topological  invariant  in  two 
dimensions,  the  relevant  action  in  this  case  has  to  be  a  curvature-squared  action  supplemented  by 
a  cosmological  constant  term.  Specifically,  part  of  the  results  in  |[23|  [22j  can  be  summarized  as 
follows.  For  a  given  lattice,  one  finds  for  the  counting  of  zero  modes 

Tetrahedron  =  4)  :  2  zero  modes 

Octahedron(lVo  =  6)  :  6  zero  modes 

Icosahedron(lVo  =  12)  :  18  zero  modes  . 

(28) 

Thus  if  the  number  of  zero  modes  for  each  regular  triangulation  of  the  sphere  is  denoted  by  Nz.m_, 
then  the  results  can  be  reexpressed  as 

Nz.m.  =  2  N0  -  6  ,  (29) 

which  agrees  with  the  expectation  that,  in  the  continuum  limit,  Nq  — >•  oo,  NzmjNo  should  approach 
the  constant  value  d  in  d  space-time  dimensions,  the  expected  number  of  local  parameters  for  a 
diffeomorphism.  Similar  estimates  were  obtained  when  looking  at  deformations  of  a  flat  lattice 
in  various  dimensions  [22J.  The  case  of  near-flat  space  is  obviously  the  simplest:  by  moving  the 
location  of  the  vertices  around  in  flat  space,  one  can  find  a  different  assignment  of  edge  lengths  that 
represents  the  same  flat  geometry.  This  then  leads  to  the  d  ■  A^-parameter  family  of  transformations 
for  the  edge  lengths  in  flat  space. 

In  general,  lattice  diffeomorphisms  correspond  to  local  deformations  of  the  edge  lengths  about 
a  vertex,  which  leave  the  local  geometry  physically  unchanged,  the  latter  being  described  by  the 

4  Instead  of  the  combinatorial  factor  n^h,  one  could  insert  a  ratio  of  volumes  Vah/Vh  (where  V%  is  the  volume  per 
hinge  nz]  and  V^h  is  the  amount  of  that  volume  in  the  simplex  a),  but  the  above  form  is  simpler. 
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values  of  local  lattice  operators  corresponding  to  local  volumes,  curvatures,  etc.  The  lesson  is  that 
the  correct  count  of  continuum  zero  modes  will,  in  general,  only  be  recovered  asymptotically  for 
large  triangulations,  where  Nq  is  significantly  larger  than  the  number  of  neighbors  to  a  point  in 
d  dimensions.  With  these  observations  in  mind,  we  can  now  turn  to  a  discussion  of  the  solution 
method  for  the  lattice  Wheeler-DeWitt  equation  in  2  +  1  dimensions. 

One  item  that  needs  to  be  discussed  at  this  point  is  the  proper  normalization  of  various  terms 
(kinetic,  cosmological  and  curvature)  appearing  in  the  lattice  equation  of  Eq.  (1211).  For  the  lattice 
gravity  action  in  d  dimensions  one  has  generally  the  following  correspondence 

f  ddx  yfg  A — y"Vg  (30) 

a 

where  Va  is  the  volume  of  a  simplex;  in  two  dimensions  it  is  simply  the  area  of  a  triangle.  The 
curvature  term  involves  deficit  angles  in  the  discrete  case, 


\  f  ddx  y/g  R  i — ■>  ^  Vh  6h  (31) 

h 

where  5h  is  the  deficit  angle  at  the  hinge  h,  and  Vh  the  associated  “volume  of  the  hinge”  |12j. 
In  four  dimensions  the  latter  is  the  area  of  a  triangle  (usually  denoted  by  Ah),  whereas  in  three 
dimensions  it  is  simply  given  by  the  length  lh  of  the  edge  labeled  by  h.  In  two  dimensions  Vh  =  1. 
In  this  work  we  will  focus  almost  exclusively  on  the  case  of  2  +  1  dimensions;  consequently  the 
relevant  formulas  will  be  Eqs.  (1301)  and  (1311)  for  dimension  d  =  2. 

The  continuum  Wheeler-DeWitt  equation  is  local,  as  can  be  seen  from  Eq.  cm ■  One  can 
integrate  the  Wheeler-DeWitt  operator  over  all  space  and  obtain 


j-  (16vrG)2 


d2x  A(g)  +  2A  J  d2Xy/g  —  J  d2Xy/gR 


^  =  0 


(32) 


with  the  super-Laplacian  on  metrics  defined  as 


Mg) 


Gij^kl  (x) 


h2 

5gij{x)5gki(x)  ' 


(33) 


In  the  discrete  case  one  has  one  local  Wheeler-DeWitt  equation  for  each  triangle  [see  Eqs.  (1201)  and 
(ED)],  which  therefore  takes  the  form 


-(16vrG)2  A(Z2)-  n^2Si  +  2\AA  l  T  =  0,  (34) 

l  iC  A  J 

where  A (Z2)  is  the  lattice  version  of  the  super-Laplacian,  and  we  have  set  for  convenience  k  =  2  nah- 
As  we  shall  see  below,  for  a  lattice  of  fixed  coordination  number,  k  is  a  constant  and  does  not  depend 
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on  the  location  on  the  lattice.  In  the  above  expression  A (l2)  is  a  discretized  form  of  the  covariant 
super-Laplacian,  acting  locally  on  the  space  of  s  =  l2  variables.  From  Eqs.  (1261)  and  (1341)  one  has 
explicitly 


A  (l2)  =  -4Aa 


dl 


d 2  d 2 

+ 


(35) 


v  da  db  db  dc  dc  da  / 

Note  that  the  curvature  term  involves  three  deficit  angles  Si,  associated  with  the  three  vertices  of  a 
triangle.  Now,  Eq.  (1M1)  applies  to  a  single  given  triangle,  with  one  equation  to  be  satisfied  at  each 
triangle  on  the  lattice.  One  can  also  construct  the  total  Hamiltonian  by  simply  summing  over  all 
triangles,  which  leads  to 


|-(16vrG)2  Y,  A(*2)  +  2A^Ha  -  «  Y  E  *  )  *  =  0  • 
l  A  A  A  icA  ) 


(36) 


Summing  over  all  triangles  (A)  is  different  from  summing  over  all  lattice  sites  (i),  and  the  above 
equation  is  equivalent  to 


j  -(16tt G)2  Y  A(^2)  +  2\Yaa  ~  Kq  Y5i  \  *  =  ° 

l  A  A  i  ) 


(37) 


where  q  is  the  lattice  coordination  number  and  is  determined  by  how  the  lattice  is  put  together 
(which  vertices  are  neighbors  to  each  other,  or,  equivalently,  by  the  so-called  incidence  matrix). 
Here,  q  is  the  number  of  neighboring  simplexes  that  share  a  given  hinge  (vertex).  For  a  flat 
triangular  lattice  q  =  6,  whereas  for  a  tetrahedron,  octahedron,  and  icosahedron,  one  has  q  =  3,  4,  5, 
respectively.  For  proper  normalization  in  Eq.  (1361)  one  requires 


d2Xy/g 



(38) 


as  well  as 


\  J  d2x  y/g  R  i — ■>  Y  Si 


This  last  correspondence  allows  one  to  fix  the  overall  normalization  of  the  curvature  term 

2 

k  =  2  na  h  —  —  , 

q 

which  then  determines  the  relative  weight  of  the  local  volume  and  curvature  terms. 


(39) 


(40) 


5  Choice  of  Coupling  Constants 

As  in  the  Euclidean  lattice  theory  of  gravity,  we  will  find  it  convenient  here  to  factor  out  an  overall 
irrelevant  length  scale  from  the  problem  and  set  the  (unsealed)  cosmological  constant  equal  to  one 
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as  was  done  in  P23-  Indeed  recall  that  the  Euclidean  path  integral  weight  always  contains  a  factor 
P(V)  oc  exp(— AoV),  where  V  =  f  yfg  is  the  total  volume  on  the  lattice,  and  Ao  is  the  unsealed 
cosmological  constant.  The  choice  Ao  =  1  then  fixes  this  overall  scale  once  and  for  all.  Since 
Ao  =  2A/167rCr,  one  then  has  A  =  8irG  in  this  system  of  units.  In  the  following  we  will  also  find  it 
convenient  to  introduce  a  scaled  coupling  A  defined  as 

A  =  I  (isb)2  <41> 

so  that  for  Ao  =  1  (in  units  of  the  UV  cutoff  or,  equivalently,  in  units  of  the  fundamental  lattice 
spacing)  one  has  A  =  1/647 tG.  One  can  now  rewrite  the  Wheeler-DeWitt  equation  so  that  the 
kinetic  term  (the  term  involving  the  Laplacian)  has  a  unit  coefficient  and  write  Eq.  (1141)  as 

(—  A  +  2A  9  - 1—^^r}^  =  0.  (42) 

I  (16ttG)2V  (16ttG)2V  J  V  ' 

Note  that  in  the  extreme  strong  coupling  limit  ( G  — >•  oo)  the  kinetic  term  is  the  dominant  one, 

followed  by  the  volume  (cosmological  constant)  term  (using  the  facts  about  A  given  above)  and, 

finally,  by  the  curvature  term.  Consequently,  at  least  in  a  first  approximation,  the  curvature  R 

term  can  be  neglected  compared  to  the  other  two  terms  in  this  limit. 

Two  further  notational  simplifications  will  be  done  in  the  following.  The  first  one  is  introduced 
in  order  to  avoid  lots  of  factors  of  167T  in  many  of  the  subsequent  formulas.  Consequently,  from 
now  on  we  shall  write  G  as  a  shorthand  for  167T  G, 


167t  G 


G 


(43) 


In  this  notation  one  then  has  A  =  G/2  and  A  =  1/4G.  The  above  notational  choices  then  lead  to  a 
much  more  streamlined  representation  of  the  Wheeler-DeWitt  equation, 

{~A  +  ~  h ^R}  ^  =  0-  (44) 

A  second  notational  choice  will  be  dictated  later  on  by  the  structure  of  the  wave  function  solutions, 
which  will  commonly  involve  factors  of  y/G.  For  this  reason  we  will  now  define  the  new  coupling  g 
as 

9  =  VG,  (45) 


so  that  A  =  4 / g2  (the  latter  g  should  not  be  confused  with  the  square  root  of  the  determinant  of 
the  metric). 

Later  on  it  will  be  convenient  to  define  a  parameter  (5  for  the  triangulations  of  the  sphere, 
defined  as 


P  = 


2  vr 

V\G2  ' 


(46) 
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Factors  of  2n  arise  here  because  we  are  looking  at  various  triangulations  of  the  two-sphere.  More 
generally,  for  a  two-dimensional  closed  manifold  with  arbitrary  topology,  one  has  by  the  Gauss- 
Bonnet  theorem 

J  d2x  yfg  R  =  4ttx  (47) 

with  x  as  the  Euler  characteristic  of  the  manifold.  The  latter  is  related  to  the  genus  g  (the  number 
of  handles)  via  %  =  2—2 g  (note  that  for  a  discrete  manifold  in  two  dimensions  one  has  the  equivalent 
form  due  to  Euler  x  =  -/Vo  —  N\  +  IV2,  where  Nt  denotes  the  number  of  simplices  of  dimension  i). 
Thus  for  a  general  two-dimensional  manifold  we  will  define 


P  = 


7T  X 

V\g2  ' 


Equivalently,  using 


V\G2  =  — =  ■  G2  =  ±G3/2 
2v  G 


and  then  making  use  of  the  coupling  g ,  one  has  simply 


p  =  % 

9 3 


for  the  sphere,  and  in  the  more  general  case 


<3  =  ^ 

9 3 


(48) 


(49) 


(50) 


(51) 


6  Outline  of  the  General  Method  of  Solution 

It  should  be  clear  from  the  previous  discussion  that  in  the  strong  coupling  limit  (large  G )  one  can,  at 
least  at  first,  neglect  the  curvature  term,  which  can  then  be  included  at  a  later  stage.  This  simplifies 
the  problem  quite  a  bit,  as  it  is  the  curvature  term  that  introduces  complicated  interactions  between 
neighboring  simplices  (this  is  evident  from  the  lattice  Wheeler-DeWitt  equation  of  Eq.  (|21l).  where 
the  deficit  angles  enter  the  curvature  term  only). 

The  general  procedure  for  finding  a  solution  will  be  as  follows.  First  a  solution  will  be  found  for 
equilateral  edge  lengths  s.  Later  this  solution  will  be  extended  to  determine  whether  it  is  consistent 
to  higher  order  in  the  weak  field  expansion.  Consequently  we  shall  write  for  the  squared  edge 
lengths 

tfj  =  s  (1  +  chij)  i  (52) 


16 


with  e  a  small  expansion  parameter.  Therefore,  for  example,  in  Eq.  (1351)  one  has  a  =  s(l  +  eha), 
b  =  s(l  +  ehb)  and  c  =  s(l  +  ehc).  The  resulting  solution  for  the  wave  function  will  then  be 
given  by  a  suitable  power  series  in  the  h  variables.  Nevertheless,  in  some  rare  cases  (such  as  the 
single-triangle  case  described  below  or  the  single  tetrahedron  in  3  +  1  dimensions  HI),  one  is  lucky 
enough  to  find  immediately  an  exact  solution,  without  having  to  rely  in  any  way  on  the  weak  field 
expansion. 

To  lowest  order  in  h,  a  solution  to  the  Wheeler-DeWitt  equation  is  readily  found  using  the 
standard  power  series  (or  Frobenius)  method,  appropriate  for  the  study  of  quantum  mechanical 
wave  equations.  In  this  method  one  first  obtains  the  correct  asymptotic  behavior  of  the  solution 
for  small  and  large  arguments  and  later  constructs  a  full  solution  by  writing  the  remainder  as  a 
power  series  or  polynomial  in  the  relevant  variable.  Of  some  importance  in  the  following  is  the 
correct  determination  of  the  wave  functional  for  small  and  large  areas  (small  and  large  s ),  and 
to  what  extent  the  resulting  wave  function  can  be  expressed  in  terms  of  invariants  such  as  areas 
and  curvatures,  or  powers  thereof. 

In  the  following  we  will  see  that  the  natural  variable  for  displaying  results  is  the  scaled  total 
area  x,  defined  as 

x  =  2 VlAtot  =  2VA  •  (53) 

A 

We  will  look  at  a  variety  of  two-dimensional  lattices,  including  the  regular  triangulations  of  the  two- 
sphere  given  by  the  tetrahedron,  octahedron  and  icosahedron,  as  well  as  the  case  of  a  triangulated 
torus  with  coordination  number  six.  In  the  equilateral  case  the  natural  variable  for  displaying  the 
results  is  then 

x  =  2  y/XAtot  =  2NaV\Aa.  (54) 

Later  on  we  will  be  interested  in  taking  the  infinite  volume  limit,  defined  in  the  usual  way  as 

oo  , 
oo  , 

const.  .  (55) 

It  follows  that  this  last  ratio  can  be  used  to  define  a  fundamental  lattice  spacing  Iq,  for  example 
via  Atot/NA  =  Aa  =  \/3Zo/4. 

The  full  solution  of  the  quantum  mechanical  problem  will,  in  general,  require  that  the  wave 
functions  be  properly  normalized,  as  in  Eq.  (1131).  This  will  introduce  at  some  stage  wave  function 
normalization  factors  A f  and  Af,  which  will  be  fixed  by  the  standard  rules  of  quantum  mechanics. 


Na  -> 

A-tot  —■ ► 
Atot 

Na 
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(56) 


If  the  wave  function  depends  on  the  total  area  only,  then  the  relevant  requirement  becomes 

r  oo  i  poo 

/  dAtot\^{Atot)\2  =  - —  /  dx  |  'F(x)  | 2  =  1. 

J o  2  V  A  Jo 

As  in  nonrelativistic  quantum  mechanics,  two  solutions  are  expected,  only  one  of  which  will  be 
regular  as  the  origin  and  thus  satisfy  the  wave  function  normalizability  requirement. 

At  this  point  it  will  be  necessary  to  discuss  each  lattice  separately  in  some  detail.  For  each 
lattice  geometry,  we  will  break  down  the  presentation  into  four  separate  items: 

(a)  Equilateral  Case  in  the  Strong  Coupling  Limit  (e  =  0).  This  subsection  will  find  a 
solution  in  the  extreme  strong  coupling  limit  (large  G ),  without  curvature  term  in  the  Wheeler- 
DeWitt  equation.  The  solution  will  not  rely  on  the  weak  field  expansion,  and  the  results  will  be 
exact  to  zeroth  order  in  the  weak  field  expansion  of  Eq.  (1521) .  In  this  case  the  simplices  are  all  taken 
to  be  equilateral,  and  the  lattice  edge  lengths  fluctuate  together. 

(b)  Equilateral  Case  with  Curvature  Term  (e  =  0).  Next,  the  curvature  term  is  included. 
The  solution  again  will  not  rely  on  the  weak  field  expansion,  and  all  the  triangles  will  be  taken  to 
be  equilateral.  The  resulting  solution  will,  therefore,  be  valid  again  (and  exact)  to  zeroth  order  in 
the  e  expansion  parameter  of  Eq.  (1521) . 

(c)  Large  Area  in  the  Strong  Coupling  Limit  (e  7^  0).  In  this  case  we  will  look  at  nonzero 
local  fluctuations  in  Eq.  (1521).  The  method  of  solution  will  now  rely  on  the  weak  field  expansion 
for  large  areas  (large  s),  but  nevertheless  it  will  turn  out  that  an  exact  solution  can  be  found  in 
this  case.  The  resulting  answer  will  be  found  to  be  correct  to  arbitrarily  large  order  0(en),  with  n 
a  positive  integer. 

(d)  Small  Area  in  the  Strong  Coupling  Limit  (e  7^  0).  Finally  we  will  look  at  the  case  of 
nonzero  fluctuations  [e  7^  0  in  Eq.  (|52D]  in  the  limit  of  small  areas  (small  s).  In  this  limit  we  will 
find  that,  in  general,  the  solution  can  be  written  entirely  in  terms  of  invariants  involving  total  areas 
and  curvatures  only  up  to  order  0(e)  or  0(e2),  depending  on  whether  a  further  symmetrization  of 
the  problem  is  performed  or  not. 

If  the  reader  is  not  interested  in  the  details  for  each  lattice,  he  can  skip  the  next  few  subsections 
and  go  directly  to  the  summary  presented  in  Sec.  (16.61). 


6.1  Single  Triangle  Case 

From  Eq.  (1341)  the  Wheeler-DeWitt  equation  for  a  single  triangle  reads 

{  (16ttG)2  4Aa  +  Wd~c  +  +  2XAa)  *(a’6’c)  =  °>  (57) 
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where  a,  6,  c  are  the  three  squared  edge  lengths  for  the  given  triangle,  and  Aa  is  the  area  of  the 
same  triangle.  Note  that  for  a  single  triangle  there  can  be  no  curvature  term.  Equivalently  one 
needs  to  solve 


d 2 


+ 


8 2 


+ 


d2 


da  db  db  dc  dc  da 


+  A  >  'h(  o,  6,  c )  =  0. 


If  one  sets 


^[a,b,c]  =  $[Aa]  , 

then  one  finds  the  following  equivalent  differential  equation 

.  d2<f>  d®  r  ,  ^ 

Aa  ,  ,n  +  2  — — —  +  16  A  Aa  $  =  0  . 

dA\  dAA 


(58) 

(59) 

(60) 


For  a  single  triangle  the  total  area  equals  the  area  of  the  single  triangle,  Atot  =  A  a-  Here  it  will 
be  convenient  to  define 

x  =  4VA  Atot  =  4\/X  Aa  (61) 

so  that  the  solution  will  be  function  of  this  variable  only.  Note  that  in  this  case,  and  in  this  case 
only,  we  will  deviate  from  the  general  definition  of  the  variable  x  given  in  Eq.  (1531).  One  can  then 
write  the  solution  to  Eq.  ()60l)  in  the  form 

Jn(x) 


4f(x)  =  AT- 


with 


so  that 


n  = 


4f(x)  =  A f  ■ 


J 


1/2 


(^4\/xAtot 


^4  VXAtot'j 

The  wave  function  normalization  constant  is  given  here  by 

Af  =  2  A1//4  . 


\  1/2 


(62) 


(63) 


(64) 


(65) 


Note  that  the  above  solution  is  exact,  and  did  not  require  in  any  way  the  weak  field  expansion. 
Two  alternate  forms  of  the  wave  function  are 


4 '(Am)  =  N 


sin  (4\/A Atot 

2  -^2tt  y/xA  tot 

AT  \J~^  exp  (-4/  V^Atot)  i-Fi  ^1, 2, 8/  V\Ato^j 


(66) 
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Here  \Fi(a,  b,  z)  is  the  confluent  hyper  geometric  functions  of  the  first  kind.  The  usefulness  of  the 
latter  representation  will  become  clearer  later,  when  other  lattices  are  considered  and  the  curvature 
term  is  included.  Expanding  the  solution  for  small  area  one  obtains 


which  shows  that  it  is  indeed  nonsingular  and,  thus,  normalizable. 


In  the  limit  of  large  areas,  a  solution  to  the  original  differential  equation  is  given  either  by  the 
asymptotic  behavior  of  the  above  Bessel  (here  sine)  function  (J),  the  same  limiting  behavior  for 
the  corresponding  Bessel  function  Y,  or  by  the  two  corresponding  Hankel  functions  (H). 


Nevertheless  among  those  four  solutions,  only  one  is  regular  and,  therefore,  physically  acceptable. 


The  calculation  for  a  single  triangle  can  be  regarded  as  a  useful  starting  point,  and  a  basic 


stepping  stone,  for  the  strong  coupling  expansion  in  1/G.  It  shows  the  physical  characteristics 


of  the  wave  function  solution  deep  in  the  strong  coupling  regime:  for  G  — >•  oo  the  coupling  term 
between  different  simplices,  which  is  caused  mainly  by  the  curvature  term,  disappears  and  one  ends 


up  with  a  completely  decoupled  problem,  where  the  edge  lengths  in  nonadjacent  simplices  fluctuate 
independently. 

6.2  Tetrahedron 

In  the  case  of  the  tetrahedron  one  has  4  triangles,  6  edges,  and  4  vertices,  and  3  neighboring 
triangles  for  each  vertex.  Let  us  discuss  again,  here,  the  various  cases  individually. 

(a)  Equilateral  Case  in  the  Strong  Coupling  Limit  (e  =  0) 

We  first  look  at  the  case  when  e  =  0  in  Eq.  (15211 ,  deep  in  the  strong  coupling  region  and  without 
the  curvature  term. 

Following  Eq.  (l53j)  we  define  the  scaled  area  variable  as 


(69) 


and  the  solution  will  be  found  later  to  be  a  function  of  this  variable  only.  For  equilateral  triangles 
the  wave  function  T  needs  to  satisfy 


9 

T"  +  -  T'  +  T  =  0  . 


(70) 


x 
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The  correct  solution  can  be  written  in  the  form 


T(x)  =  N 


J n  (  X j 


(71) 


with 


so  that 


n  2 


T(x)  =  A f 


■/ 


1/2 


(2VA 


A 


tot 


^2\/X  Atot 

The  wave  function  normalization  constant  is  given  by 

AT  =  V2X^  . 

Below  are  two  equivalent  forms  of  the  same  wave  function 

iin  ^2\/X  At, 


1/2 


(72) 


(73) 


(74) 


*(Atot)  =  A f 


V27T  \/\A 


tot 


=  AT  ]J  ^  exp  (^-  2  i  \/X  Atot  J  \F\  (1,2,4/  v  a  j± tot 
for  the  equilateral  case.  In  the  limit  of  small  area  one  obtains 


\/ A  At, 


T  =  M  \t  - 

IT 


r2  4 

1  _  ^  +  _E_  +  0  fx 6) 

6  120  y  J 


(75) 


(76) 


which  again  confirms  that  the  wave  function  is  regular  at  the  origin.  Since  one  is  solving  a  second- 
order  linear  differential  equation,  one  expects  two  solutions;  here,  one  is  singular  and  the  other 
one  is  not,  as  is  often  the  case  in  quantum  mechanics.  For  the  geometry  of  the  tetrahedron,  one 
solution  can  be  written  in  terms  of  Bessel  functions  of  the  first  kind  ( J) 

Ji/2(x)  [2  sinx 


y/x 


IT  X 


(77) 


The  Bessel  function  of  the  second  kind  (Y)  also  satisfies  the  same  differential  equation,  but  since 

Yi/2(x)  p2  cosx 


TV  X 


(78) 


this  second  solution  is  not  normalizable,  it  is  therefore  discarded  on  physical  grounds.  We  shall  see 
below  that  the  same  behavior  at  small  x  holds  also  for  the  nonzero  curvature  term.  Note  that  both 
of  the  above  solutions  are  real. 


5  There  are  also  linear  combinations  of  Bessel  functions  which  give  complex  Hankel  ( H )  functions.  These  satisfy 
the  Wheeler-DeWitt  equation  as  well;  however,  they  are  not  physically  acceptable  since  both  are  singular  at  the 
origin. 
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(b)  Equilateral  Case  with  Curvature  Term  (e  =  0) 

Next  we  include  the  effects  of  the  curvature  term.  To  zeroth  order  in  weak  field  expansions, 
when  all  edges  fluctuate  in  unison,  one  now  needs  to  solve  the  ordinary  differential  equation 

2  2  8 

&"  +  -&' - -  'k  +  =  0  ,  (79) 

x  x 

with  /?  =  27t /\f\G2  as  in  Eq.  (1461).  Since  the  deficit  angle  <5  =  n  at  each  vertex,  the  curvature 
contribution  for  each  triangle  is  k  ■  n  ■  3.  In  this  case  one  has,  therefore, 

1 


Ktetra  —  2 


(80) 


and,  therefore,  the  solution  is  given  by 
\k  ~  exp 

=  exp 


2 i  a/X Atot^j 

1 F1  | 

(-2  iV\Ato^j 

1  Fx  | 

Mot 


l  —  l 


2  4?- 
/—  ,  Zi,  ‘rk  t 

G2  VA 


i  VAA 


tot 


(81) 


in  the  equilateral  case,  up  to  an  overall  normalization  factor.  Note  that  in  this  case  one  had  to 
include  a  factor  of  Atot  /  (4Aa  )  (which  in  the  tetrahedron  case  equals  one)  in  the  imaginary  part  of 
the  first  argument  of  \  F\ . 

(c)  Large  Area  in  the  Strong  Coupling  Limit  (e  ^  0) 

Next  we  look  at  the  case  e  7^  0  in  Eq.  (1521).  In  the  limit  of  large  areas  one  finds  that  the  two 
independent  solutions  reduce  to 

exp  (±2i  a/X  Atot 


\[/  r^j 

X  — >  OO 


(82) 


exp  (zb  i  x) 

to  all  orders  in  e.  To  show  this,  one  sets  'k  =  eaAtot ,  where  Atot  is  a  sum  of  the  four  triangle  areas 
that  make  up  the  tetrahedron,  and  then  expands  the  edge  lengths  in  the  usual  way  according  to 
Eq.  (1521).  by  setting  a  =  s(l  +  eha),  etc.  Here  we  are  interested  specifically  in  the  limit  when  s  is 
large  and  e  is  small.  One  then  finds  that  the  rhs  of  the  lattice  Wheeler-DeWitt  equation  is  given 
to  Q(en)  by 


t  y/3s 


1 


a 


ct  +  4  A  )  ens 


4  2n  y/3n  n\ 

One  concludes  that  in  this  limit  it  is  sufficient  to  have 

a2  +  4  A  =  0, 

or  a  =  zb  2  i  VT,  to  obtain  an  exact  solution  in  the  limit  n 
limit  the  two  independent  wave  function  solutions  in  Eq. 
single-triangle  contributions. 


(£■ 


h)  + 


(83) 


(84) 

00.  Note  that  in  the  strong  coupling 
completely  factorize  as  a  product  of 
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(d)  Small  Area  in  the  Strong  Coupling  Limit  (e  7^  0) 

In  the  limit  of  small  area,  we  have  shown  before  that  the  solution  reduces  to  a  constant  in  the 
equilateral  case  [O  (e0)]  for  small  x  or  small  areas.  Beyond  the  equilateral  case  one  can  write  a 
general  ansatz  for  the  wave  function  in  terms  of  geometric  invariants 


and  then  expand  the  solution  in  e  for  small  s.  To  zeroth  order  in  e  we  had  the  solution  T  ~  Jn{x)/xn 
with  x  =  2  \/xAtot  and  n  =  1/2.  This  gives  in  Eq.  (1851)  70  =  0,  72  =  —  |  A  and  74  =  ^  A2.  To 
linear  order  [O  (e)]  one  finds,  though,  that  terms  appear  which  cannot  be  expressed  in  the  form 
of  Eq.  (1851).  But  one  also  finds  that,  while  these  terms  are  nonzero  if  one  uses  the  Hamiltonian 
density  (the  Hamiltonian  contribution  from  just  a  single  triangle),  if  one  uses  the  sum  of  such 
triangle  Hamiltonians,  then  the  resulting  solution  is  symmetrized,  and  the  corrections  to  Eq.  (1851) 
are  found  to  be  of  order  O  (e2) .  Then  the  wave  function  for  small  area  is  of  the  form 


\J/  rsj 


1  —  -  \  Az 

1  3  A  tot 


+ 


15 


A2  A4 


tot 


+  ... 


(86) 


up  to  terms  0(e2). 


6.3  Octahedron 

The  discussion  of  the  octahedron  proceeds  in  a  way  that  is  similar  to  what  was  done  before  for 
the  tetrahedron.  In  the  case  of  the  octahedron  one  has  8  triangles,  12  edges  and  6  vertices,  with  4 
neighboring  triangles  per  vertex.  Again  we  will  now  discuss  the  various  cases  individually. 

(a)  Equilateral  Case  in  the  Strong  Coupling  Limit  (e  =  0) 

Again  we  look  first  at  the  case  e  =  0  in  Eq.  (f52l) ,  deep  in  the  strong  coupling  region  and  without 
the  curvature  term.  Following  Eq.  (153|)  we  define  the  scaled  area  variable  as 

x  =  2  v'X  Atot  =  8x2  VX  Aa  (87) 

and  it  is  found  that  the  solution  is  a  function  of  this  variable  only.  For  equilateral  triangles  the 
wave  function  T  needs  to  satisfy 

4 

T  H - T'  +  T  =  0  .  (88) 

x 

The  correct  solution  can  be  written  in  the  form 

^  (89) 


23 


with 


n  = 


T(x)  =  M 


so  that 

J3/2  ^2  \/xAtotSj 

(2Vl4M)3/2 

The  wave  function  normalization  factor  is  given  by 

A f  =  \/l5  A1/4  . 

Equivalent  forms  of  the  above  wave  function  are 


(90) 


(91) 


(92) 


^(Atot)  =  N  23/2^5)  exp(-2z>/I  ■A-tot  )  iF,  (2,4,4z  a/a  ■A-tot 


=  A F 


cos  2 


VI 4 


2  \/27r  A  A? 


tot )  sin  (  2 

+ 


VI  Atot 


tot 


4  \/2tt  A3/2  4? 


tot 


(93) 


These  can  be  expanded  for  small  Atot  or  small  x  to  give 


T  =  A A 


3  \4r 


2  4 

i  -  —  +  —  +  0(x6 
10  280  v 


(94) 


We  note  here  again  that  both  Bessel  functions  of  the  first  (J)  and  second  (Y)  kind,  in  principle, 
give  solutions  for  this  case,  as  well  as  the  two  corresponding  Hankel  ( H )  functions.  Nevertheless, 
only  the  solution  associated  with  the  Bessel  J  function  is  regular  near  the  origin. 

(b)  Equilateral  Case  with  Curvature  Term  (e  =  0) 

Next,  we  include  the  effects  of  the  curvature  term.  Since  here  the  deficit  angle  5  =  27r/3  at 
each  vertex,  the  curvature  contribution  for  each  equilateral  triangle  is  k  ■  V  ■  3  =  27r  k.  For  the 
octahedron  one  has  in  Eq.  (l40l) 


^octa  —  2  • 


(95) 


With  the  curvature  term  one  finds 


T(4tot)  ~  exp  (-2  iV\Ato^j  1 F  |  2  - 


.  47T  Kocta 


=  exp 


(^— 2  i  Atot'j  1  F\  ^2  —  % 


VI G2 


4,  4  i 


VI A 


tot 


2  vr 


VI  G‘ 


4  4 i 


VI A 


tot 


(96) 


Note  that  in  this  case  one  had  to  include  a  factor  Atot/ (4:A/\),  which  in  the  octahedron  case  equals 
two. 

(c)  Large  Area  in  the  Strong  Coupling  Limit  (e  4  0) 
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In  the  limit  of  large  areas  the  two  independent  solutions  reduce  to 

T  ~  exp(±?'x)  ~  exp  (±  2  i  V~X  Atot]  (97) 

x  — y  oo  V  / 

to  all  orders  in  e.  In  other  words,  to  O  (en)  with  n  — >•  oo,  as  for  the  tetrahedron  case.  Note  also 
that  in  the  strong  coupling  limit  the  two  independent  wave  function  solutions  again  completely 
factorize  as  a  product  of  single-triangle  contributions. 

(d)  Small  Area  in  the  Strong  Coupling  Limit  (e  7^  0) 

In  the  limit  of  small  area,  the  solution  approaches  a  constant  in  the  equilateral  case.  To  go 
beyond  the  equilateral  case,  one  can  write  again  a  general  ansatz  for  the  wave  function,  written  in 
terms  of  geometric  invariants  as  in  Eq.  (1851).  Then  the  solution  can  be  expanded  in  e  for  small  s. 
To  zeroth  order  in  e,  the  solution  is  T  ~  Jn(x)/xn  with  n  =  3/2.  This  gives  in  Eq.  (1851)  70  =  0, 
72  =  —  |  A  and  74  =  ^  A2.  However,  to  linear  order  [ O  (e)]  one  finds  again  that  linear  terms  in 
h  appear  which  cannot  be  expressed  in  the  form  of  Eq.  (1851).  But  one  also  finds  that  while  these 
terms  are  nonzero  if  one  uses  the  Hamiltonian  density  (the  Hamiltonian  contribution  from  just 
a  single  triangle),  if  one  uses  the  sum  of  such  triangle  Hamiltonians  then  the  resulting  solution 
is  symmetrized,  and  the  corrections  to  Eq.  (1851)  are  found  to  be  of  order  O  (e2) .  Then  the  wave 
function  for  small  area  is  of  the  form 

*^1-§AAL  +  |A2  <*  +  ...  (98) 

up  to  terms  of  0(e). 


6.4  Icosahedron 

The  discussion  of  the  icosahedron  proceeds  in  a  way  that  is  similar  to  what  was  done  before  for  the 
other  regular  triangulations.  Here  one  has  20  triangles,  30  edges  and  12  vertices,  with  5  neighboring 
triangles  per  vertex.  Let  us  again  discuss  the  various  cases  individually. 

(a)  Equilateral  Case  in  the  Strong  Coupling  Limit  (e  =  0) 

Again  we  look  first  at  the  case  e  =  0  in  Eq.  ([52]) ,  deep  in  the  strong  coupling  region  and  without 
curvature  term.  Following  Eq.  (1531)  we  define  the  scaled  area  variable  as 

x  =  2  y/X  Atot  =  20x2^1  Aa  (99) 

and  a  solution  is  found  which  is  a  function  of  this  variable  only.  For  equilateral  triangles  the  wave 
function  dt  needs  to  satisfy 

4'"  +  —  T'  +  T  =  0  .  (100) 

x 
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A  solution  can  then  be  found  of  the  form 


'I'(x)  =  Af 


Jr i(®) 


(101) 


with 


so  that 


9 

n  =  — 
2 


^(x)  =  A f 


Jg/2  (2  \f\  Atot 


(2  A  am)9/2 

The  wave  function  normalization  factor  is  given  by 

Af  =  9V12155  A1/4  . 

Below  is  an  equivalent  form  of  the  same  solution 


(102) 


(103) 


4 '{Atot)  =  Af  29/2  * pjy  exp  (~2iV\Atot')  iE,  (5,  10,  4z\/XAtot 


For  small  area  Afoi  or  small  x,  one  obtains 

1 


T  =  Af 


X~  X  m  ,  fix 

1  “  ^  +  7717  +  °  (x  ) 


(104) 


(105) 


(106) 


29/2  r  (U)  [  22  1144 

which  shows  that  the  above  solution  is  regular  at  the  origin  and  normalizable. 

(b)  Equilateral  Case  with  Curvature  Term  (e  =  0) 

Next  we  include  again  the  effects  of  the  curvature  term.  Since  now  the  deficit  angle  5  =  7r/3  at 
each  vertex,  the  curvature  contribution  for  each  triangle  is  k  ■  ^  •  3  =  n  k.  For  the  icosahedron  one 
has  in  Eq.  (l40l) 


K-ir 


=  2 


1 


(107) 


Then  with  the  curvature  term  included  for  equilateral  triangles  one  obtains  for  equilateral  triangles 

(0(£°)] 


V(Ato t)  -  exp 


=  exp 


(-2;  VI  aJ)  1 F,  (5 -i^g£2,  10.  4iVX  A: 
(-2i\/I  Atot'j  1F1  (5  —  i  ^ 


-,  10,  4 i 


(108) 


up  to  an  overall  wave  function  normalization  constant.  Note  that  in  this  case  one  had  to  include  a 
factor  Afot/4AA,  which  in  the  dodecahedron  case  equals  five. 

(c)  Large  Area  in  the  Strong  Coupling  Limit  (e  7^  0) 
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In  the  limit  of  large  areas  the  two  independent  solutions  reduce  to 

T  ~  exp  (±  i  x)  ~  exp  ( ±  2  i  Atot\  (109) 

to  all  orders  in  the  weak  field  expansion  parameter  e,  as  for  the  tetrahedron  and  octahedron  case. 
Note  also  that  in  the  strong  coupling  limit  the  two  independent  wave  function  solutions  again 
completely  factorize  as  a  product  of  single-triangle  contributions. 

(d)  Small  Area  in  the  Strong  Coupling  Limit  (e  /  0) 

In  the  limit  of  small  area,  the  solution  approaches  a  constant  in  the  equilateral  case.  To  go 
beyond  the  equilateral  case,  one  can  write  again  a  general  ansatz  for  the  wave  function,  written  in 
terms  of  geometric  invariants  as  in  Eq.  (18511.  Then  the  solution  in  e  for  small  s.  To  zeroth  order 
in  e  the  solution  is  \k  ~  Jn(x)/xn  with  n  =  9/2.  This  gives  in  Eq.  (1851)  7o  =  0,  72  =  —  jj  A  and 
74  =  -j-|rj  A2.  But  to  linear  order  [O  (e)]  one  finds  again  that  linear  terms  in  h  appear  which  cannot 
be  expressed  in  the  form  of  Eq.  (|85l).  But  one  also  finds  that,  while  these  terms  are  nonzero  if 
one  uses  the  Hamiltonian  density  (the  Hamiltonian  contribution  from  just  a  single  triangle),  if  one 
uses  the  sum  of  such  triangle  Hamiltonians  then  the  resulting  solution  is  symmetrized,  and  the 
corrections  to  Eq.  (185[1  are  found  to  be  of  order  O  (e2) .  Then  the  wave  function  for  small  area  is  of 
the  form 

T  ~  1  -  yjA  A2ot  +  jfgA 2  Ajot  +  . . .  ,  (110) 

up  to  terms  of  0(e). 


6.5  Torus 

Finally  we  will  consider  a  regularly  triangulated  torus,  which  will  consist  here  of  an  infinite  lattice 
built  out  of  triangles,  with  each  triangle  having  12  neighboring  triangles.  The  torus  topology  is 
equivalent  to  requiring  periodic  boundary  conditions  in  the  two  spatial  directions.  Of  course,  one 
could  consider  the  same  type  of  lattice  but  with  some  other  sort  of  boundary  condition,  but  we 
shall  not  pursue  that  aspect  here. 

Due  to  the  local  structure  of  the  lattice  Wheeler-DeWitt  equation  in  Eq.  (fMl).  it  will  not  be 
necessary  to  include  in  the  wave  function  triangles  that  are  arbitrarily  far  apart.  Instead  it  will  be 
sufficient,  in  order  to  determine  the  overall  structure  of  the  solution,  to  include  only  those  triangles 
that  are  affected  in  a  nontrivial  way  by  the  interaction  terms  in  the  Wheeler-DeWitt  equation.  In 
the  present  case,  this  requires  the  consideration  of  one  given  triangle  plus  its  12  neighbors,  giving 
a  total  of  13  triangles.  Here,  we  will  also  set  as  before  x  =  2 


27 


(a)  Equilateral  Case  in  the  Strong  Coupling  Limit  (e  =  0) 

For  this  case  the  relevant  equation  and  its  solution  are  largely  in  line  with  what  was  obtained 
for  the  previous  cases.  For  equilateral  triangles  the  wave  function  *F  has  to  satisfy 


13 

’F"  d - V  +  >F  =  0  . 

2x 


The  wave  function  can  now  be  written  as 


T(x)  =  A f 


Jn(%) 


X“ 


with,  here,  (due  to  our  specific  choice  of  sublattice) 


11 


n  = 


so  that 


'F(x)  =  M 


J 


11/4 


^2\/a  At0t 


\  H/4  ’ 


(2  viAtoty 

The  wave  function  normalization  constant  is  given  in  this  case  by 


J 3or(1f)  -1/4 

V  r(f) 


For  the  above  wave  function  an  equivalent  form  is 


Expanding  the  above  solution  for  small  area  one  obtains 
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<F  =  A f 


2  4 

1  _  +  iL  +  o(x  6 

15  570  V 


(111) 


(112) 
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(114) 
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^(Atot)  =  U  2ll/4  ~  ^15^  exp  [~2iV\Ato^j  ^  y. ,4i  V\Ato^j  .  (116) 


(117) 


211/4  r  (^) 

which  shows  the  above  solution  is  indeed  regular  at  the  origin. 

(b)  Equilateral  Case  with  Curvature  Term  (e  =  0) 

In  the  case  of  the  torus,  the  curvature  term  is  zero  (x  =  0),  so  there  are  no  changes  to  the 
preceding  discussion. 

(c)  Large  Area  in  the  Strong  Coupling  Limit  (e  7^  0) 

In  the  limit  of  large  areas  the  two  independent  solutions  reduce  to 


X  — »■  OO 


exp  (d=  i  x)  ~  exp  (  ±  i 


i  2  VA  Atot'j 


(118) 


to  all  orders  in  e.  This  is  similar  to  what  was  found  earlier  for  the  other  lattices.  In  particular,  the 
two  independent  solutions  again  completely  factorize  as  a  product  of  single-triangle  contributions. 


28 


(d)  Small  Area  in  the  Strong  Coupling  Limit  (e  7^  0) 

In  the  limit  of  small  area,  the  regular  solution  approaches  a  constant  and  the  discussion,  and 
solution,  is  rather  similar  to  the  previous  cases.  Here  one  finds 

^  —  1  —  T5^Atot  +  2§5  A  2  Afot  +  ...  ,  (119) 

up  to  terms  of  0(e2). 


6.6  Summary  of  Results 


In  this  section  we  will  summarize  the  results  obtained  so  far  for  the  various  finite  lattices  considered 
(tetrahedron,  octahedron,  icosahedron,  and  regularly  triangulated  torus). 

(a)  Equilateral  Case  in  the  Strong  Coupling  Limit  (e  =  0) 

It  is  rather  remarkable  that  all  of  the  previous  cases  (except  the  trivial  case  of  a  single  triangle, 
which  has  no  curvature)  can  be  described  by  one  single  set  of  interpolating  wave  functions,  where 
the  interpolating  variable  is  simply  related  to  the  overall  lattice  size  (specifically,  the  number  of 
triangles) . 

Indeed  for  equilateral  triangles  and  in  the  absence  of  curvature,  the  wave  function  'k(x)  for  all 
previous  cases  is  a  solution  to  the  following  equation 

+  2n  + 1  T'  +  T  =  0  ,  (120) 

x 

with  parameter  n  given  by 

n  =  \  {NA-  2)  (121) 

where  NA  =  IV2  is  the  total  number  of  triangles  on  the  lattice.  Thus 


and,  consequently, 


Na  =  4  (n  +  i) 
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(123) 


Note  that  for  a  single  triangle  one  has  n  =  4  as  well,  but  the  definition  of  the  scaled  area  is  different 
in  that  case. 
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Furthermore  the  differential  equation  in  Eq.  (1 1201)  describes,  in  spherical  coordinates  and  with 
suitable  choice  of  constants,  the  radial  wave  function  for  a  free  quantum  particle  in  D  =  2n  +  2 
dimensions.  Indeed  recall  that  in  D  dimensions  the  Laplace  operator  in  spherical  coordinates  has 
the  form 

rFvD  n  -  1  1 

(124) 


.  d2F  D-  I9f  1 
A,f  -  +  +  ^  sD-' 


where  Agn-i  is  the  Laplace-Beltrami  operator  on  the  (D  —  l)-sphere.  In  our  case,  the  wave  function 
does  not,  to  this  order,  depend  on  angles  and  therefore  the  last  (angular  variable)  term  does  not 
contribute.  The  role  of  the  angles  is  played  in  our  case  by  the  h  variables,  which  to  this  order  do 
not  fluctuate. 

A  nonsingular,  normalizable  solution  to  Eq.  (112011  is  then  given  by 

Jn(%) 


'F(a:)  =  M  -  -  -  =  Me  lx  \F\  (n  +  \,2n  +  l,2ix) 


where  M  is  the  wave  function  normalization  constant 


M  =  2 


r(n  +  i)r(2?r  +  i) 

rfaT 


i\l1/2 


A1/4 


and 


M  = 
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M . 


(125) 


(126) 


(127) 


2"  T(n  +  1) 

Here  and  in  Eq.  (11251)  iFi(a,b~,  z)  denotes  the  confluent  hypergeometric  function  of  the  first  kind, 
sometimes  denoted  also  by  M(a,  b,  ;z).  In  either  form,  the  above  wave  function  is  real,  in  spite  of 
appearances.  The  general  asymptotic  behavior  of  the  solution  'F(x)  is  found  from  Eq.  (112011.  For 
small  x  one  has 

tf(x)  ~  xa  (128) 

with  index  a  =  0,  —2 n.  The  latter  solution  is  singular  and  will  be  discarded.  For  large  x  one  finds 
immediately 


T(x) 
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n+  o 
X  2 


—  exp (Fix) , 


(129) 


which  is  of  course  consistent  with  all  the  previous  results.  Indeed  the  other  possible  independent 
solution  of  Eq.  (112011  would  be 


’F(x)  ~ 


Yn(x) 


(130) 


where  Yn(x)  is  a  Bessel  function  of  the  second  kind  (or  Neumann  function).  However,  the  latter 
leads  to  a  wave  function  *F  which  is  singular  as  x  — >•  0, 
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T(x)  ~ - r(n)  2n  x 

7 r 


n  ^—2  n 


(131) 
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and  gives,  therefore,  a  solution  which  is  not  normalizable.  For  completeness  we  record  here  the 
small  x  (small  area)  behavior  of  the  normalized  wave  function  in  Eq.  (11251) 

1 


'F(x)  ~  A f  - 
v  ’  2nT(n  +  1) 

and  the  corresponding  large  x  (large  area)  behavior 

^2  1 


'F(x)  ~  A f 


71  Xn+  2 


nir 


cos  x  — 


(132) 


(133) 


both  of  which  reflect  well-known  properties  of  the  Bessel  functions  Jn(x). 

(b)  Equilateral  Case  with  Curvature  Term  (e  =  0) 

When  the  curvature  term  is  included  in  the  Wheeler-DeWitt  equation,  and  still  in  the  limit  of 
equilateral  triangles,  one  obtains  the  following  interpolating  differential  equation 


..  2n  +  1  ,  2/3 

<F  H - + - —  1F  +  1F  =  0  , 


(134) 

x  x 

which  now  describes  the  radial  wave  function  for  a  quantum  particle  in  D  =  2  n  +  2  dimensions, 
with  a  repulsive  Coulomb  potential  proportional  to  2 (3.  The  nonsingular,  normalizable  solution  is 
now  given  by 

’F(.x)  ~  e~lx  \F\  (n  +  |  —  i(3,  2n  +  1,  2i  x)  ,  (135) 

up  to  an  overall  wave  function  normalization  constant  J\f(n,/3).  The  normalization  constant  can 
be  evaluated  analytically  but  has  a  rather  unwieldy  form  and  will  not  be  recorded  here.  Note 
that  the  imaginary  part  (/3)  of  the  first  argument  in  the  confluent  hypergeometric  function  of 
Eq.  (I135D  depends  on  the  topology  but  does  not  depend  on  the  number  of  triangles.  In  view  of 
the  previous  discussion  the  parameter  n  increases  gradually  as  more  triangles  are  included  in  the 
simplicial  geometry.  For  the  regular  triangulations  of  the  sphere,  the  total  deficit  angle  (the  sum 
of  the  deficit  angles  in  a  given  simplicial  geometry)  is  always  47t,  so  even  if  one  writes  for  the  wave 
functional  ^[Atot,  5tot],  the  curvature  contribution  Ylh^h  is  a  constant  and  does  not  contribute  in 
any  significant  way.  Note  also  that,  in  spite  of  appearances,  the  above  wave  function  is  still  real  for 
nonzero  (3.  That  \F (x)  in  Eq.  (11351)  is  a  real  function  can  be  seen,  for  example,  from  its  definition 
via  the  power  series  expansion 


'F(x) 


1  + 


2/3 


■  x  — 


1  +  2  n-  4/3 2 
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(136) 


2n  +  1  4  +  12n  +  8n2  6  (3  +  lln  +  12n2  +  4n3) 

and  again  up  to  an  overall  normalization  factor  JV(n,{3). 

The  general  asymptotic  behavior  of  the  solution  ’F(x)  is  again  easily  determined  from  Eq.  (11341). 
For  small  x  one  has 

tf(s)  ~  (137) 
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with  again  a  =  0,  —2 n,  and,  therefore,  independent  of  the  curvature  contribution  involving  (3.  The 
second  solution  is  singular  and  will  be  discarded  as  before.  For  large  x  one  finds  immediately 


'F(x)  ~  — — p  exp{  ±i  (x  —  /31nx)  }  ,  (138) 

which  is  of  course  consistent  with  all  previous  results.  It  also  shows  that  the  convergence  properties 
of  the  wave  function  at  large  x  are  not  affected  by  the  j3  term.  A  second  independent  solution  to 
Eq.  (I134D  is  given  by 

T(x)  ~  e~lx  U  (n  +  \  —  i/3,  2n  +  1,  2ix)  ,  (139) 

where  U (a,  6, ;  z)  is  the  confluent  hypergeometric  function  of  the  second  kind  (sometimes  referred 
to  as  Tricomi’s  function).  This  second  solution  is  singular  at  the  origin,  leading  to  a  wave  function 
that  is  not  normalizable  and  will  not  be  considered  further  here. 

The  asymptotic  behavior  of  the  regular  solution  for  large  argument  z  (discussed  in  standard 
quantum  mechanics  textbooks  such  as  [25,  26|  and  whose  notation  we  will  follow  here)  can  be 
obtained  from  the  asymptotic  form  of  the  confluent  hyper  geometric  function  defined  originally, 
for  small  z,  by  the  series 


1Fi(a,b,z) 
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1  +  TT\  + 
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(140) 


It  is  common  procedure  to  then  write  i F\  (a,  6,  z)  =  W\(a,  b,  z)  +  1+2(0,  6,  z),  where  W\  and  W2  are 
separately  solutions  of  the  confluent  hyper  geometric  equation 

d2F  dF 

+?+<6-2kfa-“F  =  °-  (141> 

Then  an  asymptotic  expansion  for  \F\  (or  M)  is  obtained  from  the  following  relations: 


I+i(a,6,z) 
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T^z 
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(142) 


where 


otfd  a(a  +  l)/3(/5  +  1) 

w(a,P,z)  ~  1  +  — t  + - - 

z  -y  00  2:1!  zz  2! 


(143) 


with  the  irregular  (at  the  origin)  solution  given  instead  by  the  combination  G(a,  6,  z)  =  i  W 1  (a,  6,  z)— 
iW2(a,b,  z).  One  immediate  and  useful  consequence  of  the  above  result  is  that,  as  anticipated  be¬ 
fore,  the  behavior  of  the  regular  solution  close  to  the  origin  is  not  affected  by  the  presence  of  the 
/ 3  (curvature)  term.  In  other  words,  the  wave  function  solution  T(x)  in  Eq.  ([1351)  is  always  well 
behaved  for  small  areas  and,  therefore,  leads  to  a  perfectly  acceptable,  normalizable  solution. 
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Furthermore,  the  combination  and  properties  of  arguments  in  the  confluent  hypergeometric 
function  in  Eq.  (|135l)  allow  one  to  write  it  equivalently  as  a  Coulomb  wave  function  with  (Sommer- 
feld)  parameter  ry 

Q(V)  pl+1  ■e~ip1F1(l  +  l-  in,2l  +  2,2ip)  =Fl(r1,p ),  (144) 


where  E)(ry,  p)  denotes  the  regular  Coulomb  wave  function  that  arises  in  the  solution  of  the  quantum 
mechanical  three-dimensional  Coulomb  problem  in  spherical  coordinates  [251  [26] .  The  latter  is  a 
solution  of  the  radial  differential  equation 


d2Ft 
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i  _  ^  +  -0 


Fi  =  0 


(145) 


p  p* 

with  the  actual  radial  wave  function  then  given  by  Ri(r)  =  Fi(kr)/r.  After  comparing  the  above 
equation  with  Eq.  (11351)  one  then  identifies  p  =  x,  l  =  n—\  and  ry  =  [3.  Thus  l  =  N&/A  —  1,  where 
JVa  is  the  number  of  triangles  on  the  lattice.  The  proportionality  constant  Ci  in  Eq.  (11441)  is  given 
by  the  (Garnow)  parameter 

CM  =  2'e^Jr„(,i+„l+i,))l  •  (146) 


T{21  +  2) 

One  then  has  immediately,  from  Eq.  ([1351) .  an  equivalent  representation  for  the  regular  wave  func¬ 
tion  as 
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again  up  to  an  overall  wave  function  normalization  constant  Af(n,/3).  Again  we  note  here  that, 
on  the  other  hand,  the  irregular  Coulomb  wave  function  [usually  denoted  by  Gi{rj,p)\  is  singular 
for  small  r  and  will,  therefore,  not  be  considered  here.  Further  relevant  properties  of  the  Coulomb 
wave  function  can  be  found  in  [251 123  22-  281  [29], 

The  known  asymptotics  of  Coulomb  wave  function  [27,  28,  [29]  allow  one  to  derive  the  following 
result  for  the  wave  function  T  for  large  x 
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T(x)  ~  J\f 
with  (Coulomb)  phase  shift 

Also,  from  Eq.  (j!46l). 
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(Tn{P)  =  argT(n  +  \  +  i/3)  • 
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It  is  easy  to  check  that  the  above  result  correctly  reduces  to  the  asymptotic  expression  given  earlier 
for  'I'  in  Eq.  (11331)  in  the  limit  (3  =  0.  The  structure  of  the  wave  function  in  Eq.  (11481)  implies 

that  the  norm  is  still  finite  for  (3  ^  0,  since  the  convergence  properties  of  the  wave  function  are  not 

affected  by  the  curvature  term. 

(c)  Large  Area  in  the  Strong  Coupling  Limit  (e  ^  0) 

In  the  limit  of  large  areas  the  two  independent  solutions  reduce  to 

T  ~  exp(±ix)  (151) 

X  — »•  oo 

where  x  oc  At0t-  This  is  true  without  assuming  the  weak  field  expansion,  as  was  already  the  case 
before  (see,  in  particular,  the  section  discussing  the  tetrahedron  case). 

Consequently  in  the  strong  coupling  limit  the  two  wave  function  solutions  in  Eq.  (11511)  com¬ 
pletely  factorize  as  a  product  of  single-triangle  contributions, 

T  ~  exp  2i  V^Aa^J  ,  (152) 

A 

again  up  to  an  overall  normalization  constant.  The  above  result,  anticipated  in  [1],  was  the  basis  for 
the  variational  treatment  using  correlated  product  wave  functions  given  in  our  previous  work.  Note 
also,  in  view  of  the  result  of  Eq.  (|133l).  that  the  correct  solution,  satisfying  the  required  regularity 
condition  for  small  areas,  is  actually  a  linear  combination  of  the  above  factorized  solutions. 

(d)  Small  Area  in  the  Strong  Coupling  Limit  (e  ^  0) 

In  the  limit  of  small  area,  we  have  shown  before  in  all  cases  that  the  solution  reduces  to  a 
constant  in  the  equilateral  case  [O  (e0)]  for  small  x  or  small  areas.  To  linear  order  [O  (e)]  the 
general  result  is  still  that  linear  terms  in  h  appear  which  cannot  be  expressed  in  the  form  of 
Eq.  (1851).  But  one  also  finds  that,  while  these  terms  are  nonzero  if  one  uses  the  Hamiltonian 
density  (the  Hamiltonian  contribution  from  just  a  single  triangle),  if  one  uses  the  sum  of  such 
triangle  Hamiltonians  then  the  resulting  solution  is  symmetrized,  and  the  corrections  to  Eq.  (1851) 
are  found  to  be  of  order  O  (e2) .  In  other  words,  it  seems  that  some  residual  lattice  artifacts  that 
survive  at  very  short  distances  can  be  partially  removed  by  a  suitable  coarse-graining  procedure  on 
the  Hamiltonian  density. 

One  might  wonder  what  lattices  correspond  to  values  of  n  greater  that  9/2,  which  is  the  highest 
value  attained  for  a  regular  triangulation  of  the  sphere,  corresponding  to  the  icosahedron.  For 
each  of  the  three  regular  triangulations  with  Nq  sites  one  has  for  the  number  of  edges  N±  =  |fVo 
and  for  the  number  of  triangles  IV2  =  (|  —  1)Nq  +  2,  where  q  is  the  number  of  edges  meeting 
at  a  vertex  (the  local  coordination  number).  In  the  three  cases  examined  before  q  was  between 
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three  and  five,  with  six  corresponding  to  the  regularly  triangulated  torus.  Note  that  for  a  sphere 
Nq  —  IVi  +  IV2  =  2  always.  The  interpretation  of  other,  even  noninteger,  values  of  q  is  then  clear. 
Additional  triangulations  of  the  sphere  can  be  constructed  by  considering  irregular  triangulations, 
where  now  the  parameter  q  is  interpreted  as  an  average  coordination  number.  Of  course  the 
simplest  example  is  a  semiregular  lattice  with  Na  vertices  with  coordination  number  qa  and  iV& 
vertices  with  coordination  number  q\ >,  such  that  Na  +  N),  =  Nq.  Various  irregular  and  random 
lattices  were  considered  in  detail  some  time  ago  in  m,  and  we  refer  the  reader  to  this  work  for  a 
clear  exposition  of  the  properties  of  these  lattices. 

We  conclude  this  section  by  briefly  summarizing  the  key  properties  of  the  gravitational  wave 
function  given  in  Eqs.  ([1351)  and  (11471).  which  from  now  on  will  be  used  as  the  basis  for  additional 
calculations.  First  we  note  that  the  above  wave  function  is  a  function  of  the  total  area  and  total 
curvature  only  and,  as  such,  is  manifestly  diffeomorphism-invariant  and  in  accord  with  the  spatial 
diffeomorphism  constraint.  While  it  was  derived  by  looking  at  the  discrete  triangulations  of  the 
sphere,  it  contains  a  parameter  n,  related  to  the  total  number  of  triangles  on  the  lattice  by  Eq.  (11211). 
that  will  allow  us  to  go  beyond  the  case  of  a  finite  lattice  and  investigate  the  physically  meaningful, 
and  presumably  universal,  infinite  volume  limit  n  — >•  00  [see  Eq.  (1551)].  We  have  also  shown  that  the 
above  wave  function  is,  in  all  cases,  an  exact  solution  of  the  full  lattice  Wheeler-DeWitt  equation  of 
Eq.  (1211)  in  the  limit  of  large  areas,  and  to  all  orders  in  the  weak  field  expansion.  Again,  this  last  case 
is  most  relevant  for  taking  the  infinite  volume  limit,  defined  previously  in  Eq.  (1551).  Furthermore, 
the  small  area  behavior  of  the  wave  function  plays  a  crucial  role  in  uniquely  constraining,  through 
the  regularity  condition,  the  correct  choice  of  solution.  In  this  last  limit  one  also  finds  that  the 
various  individual  lattice  solutions  agree  with  the  universal  form  of  Eqs.  (11351)  and  (11471)  only 
to  a  low  order  in  the  weak  field  expansion,  which  is  expected  given  the  different  short  distance 
lattice  artifacts  of  the  regular  triangulation  solutions.  Nevertheless,  knowledge  of  their  behavior 
is  completely  adequate  for  extracting  the  most  important  physically  relevant  piece  of  information, 
namely  the  constraint  on  the  wave  function  based  on  the  stated  regularity  condition  at  small  areas, 
which  comes  down  to  a  simple  integrability  or  power  counting  argument. 


7  Average  Area 

In  this  section  we  will  look  at  a  natural  quantum  mechanical  expectation  value,  the  average  total 
physical  area  of  the  lattice  simplicial  geometry.  It  is  one  of  many  quantities  that  can  be  calcu- 
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Octahedron 


Figure  4:  Wave  Function  \F  versus  total  area  for  the  octahedron  lattice,  with  and  without  curvature 
contribution.  The  wave  function  is  shown  here  for  g  =  \[G  =  1,  a  value  chosen  here  for  illustration 
purposes.  The  relevant  expression  for  the  wave  function  is  given  in  Eq.  (11531).  We  refer  to  the  text 
for  further  details  on  how  the  wave  function  was  obtained,  and  what  its  domain  of  validity  is.  The 
wave  functions  shown  here  have  been  properly  normalized.  Note  that  with  a  nonzero  curvature 
term  the  peak  in  the  wave  function  moves  away  from  the  origin. 


lated  within  the  lattice  quantum  gravity  formalism,  and  is  clearly  both  manifestly  geometric  and 
diffeomorphism-invariant.  Here  we  will  use  the  wave  functions  given  in  Eqs.  (11351)  and  (114711 .  origi¬ 
nally  obtained  for  the  tetrahedron,  octahedron  and  icosahedron,  and  later  extended  to  any  number 
of  triangles  Na 

xV(Atot)  ~  e~l~  iFi  +  \  -  i/3,2n  +  1,2*  >  (153) 

with  n  =  j  (Na  —  2),  /3  =  Air/g 3  and  g  =  VG,  and  again  valid  up  to  an  overall  wave  function 
normalization  constant.  Due  to  the  structure  of  the  wave  function  the  resulting  probability  dis¬ 
tribution  for  the  area  is  rather  nontrivial,  having  many  peaks  associated  with  the  infinitely  many 
minima  and  maxima  of  the  hyper  geometric  function.  Clearly  the  most  interesting  limit  is  one 
where  one  considers  an  infinite  number  of  triangles,  Na  — >  oo,  which  corresponds  to  n  — >•  oo  in 
Eq.  (|153D.  In  Figs.  [4] and [5]  we  display  the  behavior  of  the  wave  function  in  Eq.  (11531).  both  with 
and  without  the  curvature  contribution  in  the  Wheeler-DeWitt  equation.  One  notices  that  when 
the  curvature  term  is  included  (/3  7^  ()),  the  peak  in  the  wave  function  shifts  away  from  the  origin. 
This  is  largely  expected,  based  on  the  contribution  from  the  repulsive  Coulomb  term  in  the  wave 
equation  of  Eq.  (11341) . 
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Figure  5:  Same  wave  function  'll  as  in  Fig.  (U  but  now  for  the  icosahedron  lattice. 


The  average  total  area  can  then  be  computed  from  the  above  wave  function,  as  the  ground  state 
expectation  value 

,  IJMM  li'tg)!2  ,,,,, 

<»l»>  /  <M<i]  |*(s)|*  ’  1  ’ 

where  g ,  here,  is  the  three-metric,  and  dg[g\  denotes  a  functional  integration  over  all  three- metrics. 
In  our  case  we  use  the  measure 


I  iM  ->  f 


dA 


tot  ? 


(155) 


(156) 


which  then  gives,  in  terms  of  the  scaled  area  variable  x, 

^  ^  /o°°  dxx-\^>(x)\2 

“  9  r*i®wia  ' 

In  the  absence  of  a  curvature  term  in  the  Wheeler-DeWitt  equation  (/3  =  0),  the  average  area  can 
easily  be  computed  analytically  in  terms  of  Bessel  function  integrals,  and  the  result  is 

7 r  (4 n  —  1)  T(4n  —  2) 


<  Atot  >  =  g  • 


(157) 


28"-5r  (n)4 

Note  that  the  average  area  diverges  as  n  — >•  which  corresponds  to  the  tetrahedron;  this  entirely 
spurious  divergence  prevents  us  from  using  the  tetrahedron  lattice  in  plotting  and  numerically 
extrapolating  the  remaining  two  lattices  (octahedron  and  icosahedron)  to  the  infinite  lattice  limit. 
For  the  octahedron  one  finds  <  Atot  >=  15  <7/%,  for  the  icosahedron  <  Atot  >=  21879g/39207r,  and 


in  the  large  n  limit  <  Atot  >=  \J2npi r  g  +  0(\/y/n). 
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One  finds  that  in  the  presence  of  a  curvature  term  {(3  ^  0)  the  resulting  integrals  are  significantly 
more  complicated.  We  have,  therefore,  resorted  to  a  number  of  tools,  which  include  an  analytic 
expansion  in  f3,  the  use  of  known  asymptotic  expansions  for  the  wave  function  at  large  arguments, 
and  an  exact  numerical  integration  of  the  resulting  integrals.  Let  us  first  discuss  here  the  expansion 
in  f3.  It  is  known  that  the  Coulomb  wave  functions  can  be  expanded  in  terms  of  spherical  Bessel 
functions  (Neumann  expansion)  [271 1251  [29] .  so  that  one  has 


oJ+l 

Fi(ri,p)  =  + 


7 r 


Y2  bk^Jk+ 


(158) 


k=i 


with  coefficients  6^(77)  given  by  a  simple  recursion  relation.  When  written  out  explicitly,  the  ex¬ 
pression  in  curly  brackets  involves 


J 


1+  2 


.  2Z  +  3  ,  .  2/  +  5  2  t 

1 0)  +  ,  ,  ,  V-J,,  3  (x)  +  ,  ;  ;  v  ■  J. 


I  +  1 


1+  2 


Z  +  1 



(x)  + 


(159) 


with  additional  terms  linear  in  rj  reappearing  at  higher  orders.  That  the  above  expansion  is  a  bit 
problematic  is  not  entirely  surprising,  given  the  modified  asymptotic  behavior  of  the  Coulomb  wave 
functions  for  r]  ^  0.  In  the  following,  in  order  to  provide  initially  some  insight  into  the  effects  of 
the  (or  (3)  term  on  the  wave  function  T,  we  will  include  the  first  correction  as  a  perturbation, 
and  drop  the  rest.  Later  on,  higher  order  corrections  can  be  included  as  additional  contributions. 
With  this  truncation,  the  Coulomb  wave  function  in  Eq.  (11441)  becomes 


pi(v,p)  =  ^rr(*  +  !)  Cl(v) 


t  /  \  2/ +  3  .. 

J,+  i^  +  riJTTJt+l(p)  +  '" 


(160) 


with  the  last  term  treated  as  a  perturbation,  giving  for  the  wave  function  itself  [see  Eq.  ([13511] 


'L(x)  ~  e  lx  1F1  (n  +  ^  —  i(3,  2n  +  1,  2ix) 

2 rfi  — j—  2 

Jn(x )  T  (3  1  Jn+l{x~)  + 

n  +  ^ 


(161) 


again  up  to  an  overall  wave  function  normalization  constant  J\f.  Note  that  if  m  Bessel  function 
terms  are  kept  in  Eq.  (11611).  beyond  the  zeroth  order,  strong  coupling,  term  involving  Jn(x),  then 
the  resulting  expansion  in  (3  contains  terms  up  to  (3m .  One  finds  to  lowest  order  ( m  =  1) 

1  T(n)  _  Al~n{n  +  1)P 

JP  ~  7 


+ 


+ 


(162) 


2r(n  +  i)T(2n  +  i)  (2n  +  l)r(n  +  l)2 

From  the  above  expressions,  the  average  area  can  then  be  computed  as  some  still  rather  complicated 
function, 

41-2»r(n-i)r(„  +  i)r(2„  +  i)2i 


,  „  |  7T  (4n  -  1)  T(4n  -  2)  ,  4(n  +  l)/3 

<  A-tot  >—  9  \  - \  a - r 


2Sn-5  p(n)4 


2n  +  1 


1  - 


n2r(n) 


+  ■■ 
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<Aa> 


Figure  6:  Average  area  of  a  single  triangle  vs.  g  =  \[G  for  the  octahedron  and  the  icosahedron 
configurations.  The  average  area  was  calculated  using  the  expression  in  Eq.  (11561).  Note  the 
qualitative  change  when  one  includes  the  curvature  term,  with  a  minimum  appearing  at  g  rs j  0(1). 


Additional  terms  can  later  be  included  in  the  Bessel  function  expansion  of  Eq.  (11581) ,  so  as  to  obtain 
more  accurate  values  for  the  averages;  this  will  be  done  later. 

Figure  [6]  shows  the  exact  value  of  the  average  area  for  a  single  triangle  <  Aa  >=<  Atot  >  /Na 
as  a  function  of  the  coupling  g,  obtained  by  doing  the  integral  in  Eq.  (I156D  numerically,  with  the 
wave  function  given  in  Eq.  (11531).  One  noteworthy  aspect  is  that  a  qualitative  change  seems  to 
occur  when  one  includes  the  curvature  term:  a  well  defined  minimum  occurs  at  g  ~  1,  which  would 
suggest  the  appearance  of  some  sort  of  phase  transition.  Doing  the  integrals  numerically  one  finds 
a  minimum  in  the  average  area  of  a  triangle  at  gc  ~  3.1  for  the  octahedron,  and  at  gc  ~  2.6  for  the 
icosahedron.  On  the  other  hand,  using  the  lowest  order  Bessel  function  expansion  of  Eq.  (I161D  for 
the  octahedron  (n  =  3/2)  one  finds  a  minimum  at  gc  =  2.683,  and  for  the  icosahedron  (n  =  9/2)  at 
gc  =  2.271.  Adding  one  more  Bessel  function  correction  term  then  gives  gc  =  3.135  and  gc  =  2.637 
for  the  two  cases,  respectively,  which  suggests  that  the  expansion  is  converging. 

The  limit  of  a  large  number  of  triangles  Na  — >•  oo  corresponds  to  taking  the  parameter  n  in 
Eq.  (|153D  to  infinity,  since  n  =  1  (Na  —  2).  From  the  lowest  order  Bessel  function  expansion  one 
obtains  the  following  analytic  expression  for  the  average  total  area 


<  Atot  > 


+ 


2(tt-2  )g 
7 r 


P  +  '  "  5 


(164) 


with  f3  =  4-7T /g3  [see  Eq.  (j50|)].  In  this  limit  the  resulting  function  of  g  has,  again,  a  well  defined 
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minimum  at 


(165) 


o  8(ir-2)V2^ 

9c  =  - 7 = - 

y/n 

or  gc  ~  2.839 /n1/6  for  large  n  with  one  Bessel  function  correction  term.  With  two  Bessel  function 
correction  terms  in  Eq.  (|161l)  one  finds  gc  —  3.276 /n1//6,  which  again  suggests  that  the  expansion 
is  slowly  converging.  Using  the  exact  wave  function  to  do  the  integrals  numerically  one  finds 
for  the  minimum  gc  ~  3.309/n1/6,  which  is  close  to  the  above  answer.  Interestingly  enough,  the 
above  result  would  suggest  that  in  the  limit  of  infinitely  many  lattice  points  the  critical  point  gc 
actually  moves  to  the  origin,  indicating  a  phase  transition  located  at  exactly  g  =  0  (G  =  0)  in  the 
infinite  volume  (n  — >•  oo)  limit  (see  further  discussion  later).  We  note  here  that  the  average  area 
for  a  single  triangle  is  obtained  by  simply  dividing  the  average  total  area  by  the  total  number  of 
triangles  N&  =  4n  +  2,  which  then  gives  in  the  same  limit  of  large  n  and  strong  coupling 

<  Aa  >  =  — t==  +  0{~)  .  (166) 

2v/2vrn  n 

Quite  generally,  the  average  of  the  area  per  site  in  the  lattice  theory  (the  spatial  volume  per  site) 
appears  to  be  well  defined  mainly  due  to  our  wave  function  normalization  choices  and,  consequently, 
can  be  explicitly  calculated  without  any  leftover  ambiguity. 

As  will  be  discussed  further  below  in  more  detail,  the  estimate  for  the  critical  point  given  in 
Eq.  (11651)  is  also  in  good  agreement  with  a  previous  variational  estimate.  In  [1]  the  quantum- 
mechanical  variational  (Rayleigh-Ritz)  method  was  used  to  find  an  approximation  for  the  ground 
state  wave  function,  using  as  variational  wave  function  a  correlated  (Jastrow-Slater)  product  of 
single-triangle  wave  functions.  There  it  was  found,  from  the  roots  of  the  equation  <  T|i7|T  >  =  0, 
that  the  variational  parameters  are  almost  purely  imaginary  for  strong  coupling  (large  G  >  Gc), 
whereas  for  weak  enough  coupling  (small  G  <  Gc )  they  become  real.  This  abrupt  change  in 
behavior  of  the  wave  function  at  Gc  then  suggested  the  presence  of  a  phase  transition.  With  the 
notation  used  in  this  paper,  the  result  of  [I]  reads  g').  ~  1/N&,  in  qualitative  agreement  with  the 
result  of  Eq.  (|165j).  in  the  sense  that  both  calculations  point  to  a  critical  point  Gc  =  0  in  the  infinite 
volume  limit. 

Let  us  now  make  some  additional  comments  which  should  help  clarify  the  interpretation  of  the 
previous  results.  It  is  well  known  that  if  there  is  some  sort  of  continuous  phase  transition  in  the 
lattice  theory,  the  latter  is  generally  associated  with  a  divergent  correlation  length  in  the  vicinity 
of  the  critical  point.  In  our  case  it  is  clear  that  at  strong  coupling  (large  g)  the  correlation  length  is 
small  (of  order  one)  in  units  of  lattice  spacing.  This  can  be  seen  from  the  fact  that  (a)  the  coupling 
term  in  the  Wheeler-DeWitt  equation  is  due  mainly  to  the  curvature  term,  which  is  small  for  large 
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g,  and  (b)  that  the  ground  state  wave  function  is  of  the  form  of  a  correlated  product  in  the  same 
limit  [see  Eq.  ([15211] .  Then  as  the  effects  of  the  curvature  term  are  included,  the  correlation  length 
starts  to  grow  due  to  the  additional  coupling  between  edge  variables.  The  previous  calculation 
would  then  suggest  that  the  point  of  divergence  is  located  at  g  =  0.  It  is,  of  course,  essential  that 
one  looks  at  the  limit  of  infinitely  many  triangles,  N&  — >•  oo,  since  no  continuous  phase  transition 
can  occur  in  a  system  with  a  finite  number  of  degrees  of  freedom. 

It  is  also  of  interest  here  to  discuss  how  the  above  (Lorentzian)  results  relate  to  what  is  known 
about  the  corresponding  Euclidean  lattice  theory  in  three  dimensions,  which  was  studied  in  some 
detail  in  [23].  There  a  phase  transition  was  found  between  two  phases,  with  the  weak  coupling  phase 
G  <  Gc  exhibiting  a  sort  of  pathological  behavior,  whereby  the  lattice  collapses  into  what  geomet¬ 
rically  could  be  described  as  a  branched  polymer.  This  is  clearly  a  nonperturbative  phenomenon 
that  cannot  be  seen  from  perturbation  theory  in  G.  In  the  Euclidean  formulation,  average  volumes 
are  obtained  as  suitable  derivatives  of  log  Ziatt  with  respect  to  the  bare  cosmological  constant  Ao, 
where  Ziatt  is  the  lattice  path  integral 


with,  in  four  dimensions,  the  action  given  by 

Iiatt  =  Ao  X]  Vh(l2)  ~  k  Y  $h(l2)Ah(l2)  (168) 

h  h 

and  h  denoting  a  hinge  [more  details  can  be  found  in  [24]].  Similarly,  the  average  curvature  can 
also  be  obtained  as  a  derivative  of  log  Z  with  respect  to  k  =  1/(87 tG).  More  importantly,  a 
nonanalyticity  in  Z,  as  induced  by  a  phase  transition,  is  expected  to  show  up  in  local  averages  as 
well.  From  the  above  expression  for  exact  sum  rules  can  be  derived  relating  various  averages 
[30] .  In  the  case  of  the  three-dimensional  Euclidean  theory  the  sum  rule  reads 

2  Ao  <  ^  Vt  >  —k  <  $hlh  >  —  Co  =  0  (169) 

T  h 

where  the  first  term  contains  a  sum  over  all  lattice  tetrahedra,  and  the  second  term  involves  a  sum 
over  all  lattice  hinges  (just  edges  in  this  case).  The  quantity  Co,  here,  is  a  constant  that  solely 
depends  on  how  the  lattice  is  put  together  (i.e.  on  the  local  coordination  number,  or  incidence 
matrix) . 

In  [24]  it  was  found  that  the  average  curvature  goes  to  zero  at  some  gc  with  a  characteristic 
universal  exponent  6, 

<  Y  Mh  >  =  -Rq  | g  -  gc\5  (170) 

h 
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and  that  the  curvature  fluctuation  diverges  in  the  same  limit.  From  the  sum  rule  in  Eq.  (11691)  one 
then  deduces  that  the  average  volume  in  the  Euclidean  theory  has  a  singularity  of  the  type 

<J2vt>=  vo  -  Vi\g-9c\5  (171) 

T 

with  the  same  exponent  5  —  0.77.  The  latter  is  related  by  standard  universality  and  scaling 
arguments  EHEHIH]  (see  [19]  for  details  specific  to  the  gravity  case)  to  the  correlation  length 
exponent  v  by  v  =  (1  +  6)/d  in  d  dimensions.  To  compare  to  the  Lorentzian  theory  discussed 
in  this  paper,  one  notes  that  the  three-dimensional  Euclidean  theory  corresponds  to  the  (2  +  1)- 
dimensional  Wheeler-DeWitt  theory,  so  that  the  average  volume  in  the  above  discussion  should  be 
taken  to  correspond  to  an  average  area  in  our  case.  |f|  To  conclude,  the  results  for  the  average  area 
suggest  the  existence  of  a  phase  transition  in  the  Lorentzian  theory  located  at  g  =  0.  In  the  next 
sections  we  will  present  a  further  test  of  this  hypothesis,  based  on  physical  observables  that  can 
establish  directly  and  unambiguously  the  location  of  the  phase  transition  point. 


8  Area  Fluctuation,  Fixed  Point  and  Critical  Exponent 


Another  quantity  that  can  be  obtained  readily  from  the  wave  function  T  is  the  fluctuation  in  the 
total  area 

XA  =  ^  {<  {Atot?  >  ~  <  Atot  >2}  .  (172) 

The  latter  is  related  to  the  fluctuations  in  the  individual  triangles  by 

XA  =  Na{<A2a>  -  <Aa  >2}  (173) 

with  the  usual  definition  of  averages,  such  as  the  one  given  in  Eq.  (|154p. 

Generally  for  a  field  <j)(x )  with  renormalized  mass  m  and  correlation  length  £  =  m—  ,  wave 
function  renormalization  constant  Z .  and  (Euclidean)  propagator 

<  >=  [-pL  -p—  ,  (174) 

J  yZlTja  Pz  +  rnz 

one  has  for  <3?  =  J  4>(x) 

<  $2  >  =  f  <  <j>(x)<j>(y)  >  =  V  f  <  (j>(x)(j>(0)  >=  V  =  V  Zf  .  (175) 

J  x,y  Jx  ^ 

6  It  should  be  noted  that  in  the  case  of  the  lattice  Wheeler-DeWitt  equation  of  Eqs.  m  and  m  and,  generally, 
in  any  lattice  Hamiltonian  continuous-time  formulation,  the  lattice  continuum  limit  along  the  time  direction  has 
already  been  taken.  This  is  due  to  the  fact  that  one  can  view  the  resulting  2  +  1  theory  as  originating  from  one  where 
there  exist  initially  two  lattices  spacings,  at  and  a.  The  first  one  is  relevant  for  the  time  direction  and  the  second  one 
for  the  spatial  directions.  In  the  present  lattice  formulation  the  limit  at  — »  0  has  already  been  taken;  the  only  limit 
left  is  a  — >  0,  which  requires  the  existence  of  an  ultraviolet  fixed  point  of  the  renormalization  group. 
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Thus  the  field  fluctuation  probes  the  propagator  at  zero  momentum,  which  in  turn  is  directly 
related  to  the  renormalized  mass  (and  thus  £)  for  the  field  in  question.  If  the  field  $  acquires  a 
nonzero  expectation  value,  the  above  result  is  modified  to 

^  {<  $2  >  -  <  $  >2}  =  =  Z  i2  ,  (176) 

V  K  J  mz 


involving  instead  the  connected  propagator.  In  the  gravity  case  the  quantity  Atot  plays  the  role  of 
d>;  if  the  fluctuation  diverges  (£  — >•  oo)  then  one  has  a  phase  transition  or  an  ultraviolet  fixed  point 
in  quantum  field  theory  language  mm- 

Without  the  curvature  term  in  the  Wheeler-DeWitt  equation  [/3  =  0  for  the  wave  function  in 
Eq.  (11611)].  the  area  fluctuation  does  not  diverge,  even  when  n  is  large  and  is  simply  proportional 
to  g2.  In  this  case  one  finds 


Xa(P  =  0)  = 
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4n  —  1  2n  —  1  7r2(4n  —  l)T(4n  —  2) 

16  ~  2T^ZT3(2^TT)r(n)'8_ 


9 2  ~ 


^—^g2  +  0(~).  (177) 

47T  n 


Note  the  spurious  singularity  for  the  special  case  of  the  tetrahedron,  n  =  1/2.  When  the  curvature 
term  is  taken  into  account  one  finds,  from  the  full  wave  function  T  in  Eq.  (11611)  and  in  the  limit  of 
large  n, 

XA  =  (l--\  g-  +  +  ■■■  (178) 

\  7T  J  4  V  Tl  7T  g 

Note  that  the  fluctuation  now  appears  to  diverge  as  g  0  (see  also  Fig.  0).  Furthermore,  XA  is 
nonanalytic  in  the  original  Newton’s  coupling  G  =  g2,  which  suggests  that  perturbation  theory  in 
G  is  useless.  A  divergence  of  the  fluctuations  as  g  -A  0  implies  that  in  this  limit  the  correlation 
length  diverges  in  lattice  units,  signaling  the  emergence  of  a  massless  excitation. 

Just  as  for  the  case  of  the  average  curvature  [Eq.  (1169|)] .  an  exact  sum  rule  can  be  derived  in 
the  (Euclidean)  lattice  path  integral  formulation,  relating  the  local  volume  fluctuations  to  the  local 
curvature  fluctuations.  In  the  three-dimensional  Euclidean  path  integral  theory  the  following  exact 
identity  holds  for  the  fluctuations  m 


4A2 


<(X>*) 


2  >  -  <  E  Vh  >2  -  k2  <  (E  <W2  >  -  <  E  Wh  >2  -  2iW  =  0  ,  (179) 

h  _  _  h  h  _ 

where  Aq  is  the  number  of  edges  on  the  lattice  (further  exact  sum  rules  can  be  derived  by  considering 
even  higher  derivatives  of  the  free  energy  lnZ^  with  respect  to  the  parameters  Ao  and  k).  Since  the 
last  equation  relates  the  fluctuation  in  the  curvature  to  fluctuations  in  the  volumes,  it  also  implies 


a  relationship  between  their  singular  (divergent)  parts 


.0 


We  noted  previously  that  in  our  Hamiltonian  formulation,  the  lattice  continuum  limit  along  the  time  direction 
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Octahedron 



Figure  7:  Area  fluctuation  \A  vs  g  =  VG  for  the  octahedron  and  icosahedron,  computed  from 
Eq.  (I172D.  Note  the  divergence  for  small  g. 

According  to  the  sum  rule  of  Eq.  ([1791)  a  divergence  in  the  curvature  fluctuation 

XR  ~  <  <£  <5/A)2  >  —  <  >2  (180) 

h  h 

for  the  three-dimensional  (Euclidean)  theory  generally  implies  a  corresponding  divergence  in  the 
volume  fluctuation 

XV  ~  <  (E  Vh?  >  ~  <  E  V*  >2  (181) 

h  h 

for  the  same  theory.  In  our  case  a  divergence  is  expected  in  2  +  1  dimensions  of  the  form 

XA  ~  \9-9c\~a  (182) 

9^9c 

with  exponent  a  =  1  —  S  =  2  —  3r,  where  5  is  the  universal  curvature  exponent  defined  previously 
in  Eq.  (11701).  and  v  the  correlation  length  exponent.  The  latter  is  defined  in  the  usual  way  [311, [32] 
through 

£  ~  1 9-9c\~v,  (183) 

9^9c 

where  ^  is  the  invariant  gravitational  correlation  length.  The  scaling  relations  among  various 

exponents  (v,  6,  a)  are  rather  immediate  consequences  of  the  scaling  assumption  for  the  singular 

has  already  been  taken.  This  results  in  two  lattice  spacings,  one  for  the  time  and  one  for  the  space  directions,  denoted 
here  respectively  by  at  and  a,  with  the  first  lattice  spacing  already  sent  to  zero.  It  is  then  relatively  straightforward 
to  relate  volumes  between  the  two  formulations,  such  as  V  —  atA.  Relating  curvatures  (for  example,  2R  in  the  2  +  1 
theory  vs  the  Ricci  scalar  R  in  the  original  three-dimensional  theory)  in  the  two  formulations  is  obviously  less  easy, 
due  to  the  presence  of  derivatives  along  the  time  direction. 
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part  of  the  free  energy,  Fsing  ~  £  d  in  the  vicinity  of  a  critical  point  (for  more  detailed  discussion  see, 
for  example,  usual  si).  The  preceding  argument  then  implies,  via  scaling,  that  a  determination 
of  a  provides  a  direct  estimate  for  the  correlation  length  exponent  v  defined  in  Eq.  (11831).  Note 
that  based  on  the  results  so  far  one  would  be  inclined  to  conclude  that  for  2  +  1  gravity  the  critical 
point  c/c  — ^  0  as  ti  — ^  oo .  Equation  (I182D  can  then  be  rewritten  either  as 

XA  ~  e/v  (184) 

g^9c 

or,  in  a  finite  volume  with  linear  lattice  dimensions  L  ~  N^d  ~  ~  y/n  (since  IVa  =  4n  +  2), 

as 

XA  ~  La/U  ~  jr1^"3/2  ,  (185) 

9^>9c 

since,  for  a  very  large  box  and  g  very  close  to  the  critical  point  gc,  the  correlation  length  saturates 
to  its  maximum  value  £  ~  L.  Hence  the  volume-  or  n-dependence  of  x  provides  a  clear  and  direct 
way  to  estimate  the  critical  correlation  length  exponent  v  defined  in  Eq.  (11831). 


9  Results  for  Arbitrary  Euler  Characteristic  y 


The  results  of  the  previous  sections  refer  to  regular  triangulations  of  the  sphere  (x  =  2)  and  the 
torus  (x  =  0)  in  2  +  1  dimensions.  It  would  seem  that  one  has  enough  information  at  this  point 
to  reconstruct  the  same  type  of  answers  for  arbitrary  In  particular  one  has  for  the  parameter  (3 
[see  Eqs.  (1481)  and  (T5TT)] 

9-77-V 

(186) 



5° 


relevant  for  the  wave  functions  in  Eqs.  (11351)  or  (11471).  For  the  average  total  area  one  then  finds, 
using  the  wave  function  expansion  in  Eq.  (116111 . 


+  A tot  O —  9  ' 


21~2,T  (n  -  ±)  T  (2n  +  \)  ^  +  1^7r  X 


1  - 


41-2nr(n_i)r(„.+  i)r(2n+i)2 

n2r(n )6 


T(n); 


g3  (2n  +  1) 


+ 
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In  the  large  n  limit  one  obtains  for  the  average  area  of  a  single  triangle 


<  Aa  >  = 
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2  y/2- 


7 r  n 


i-a+0+ 


16n 


rr 
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(n  -  2)  x 

g  '2  n 


4n(7r  —  2) 


+  ... 


and  for  the  average  total  area 


„  „  .  ^  ,  4(tt-2)x 

<  Aot  >  ~  — y  g  +  - ^ -  + 


7 r 


(187) 


(188) 


(189) 
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For  the  area  fluctuation  defined  in  Eq.  (117311  one  finds  in  the  same  large  n  limit 

J2 


IT 


xa  =  (1-z)^r  +  c’(r)  +  ^-7T)\l-rz~  + 


n 


2  X 


mr  g 


(190) 


Again  note  that  the  fluctuation  appears  to  diverge  as  g  -A  0,  which  implies  that  this  is  the  more 
interesting  limit,  so  from  now  on  we  will  focus  specifically  on  this  limit.  It  is  clear  from  the  analytic 
expression  for  <  Atot  >  in  Eqs.  (118711  or  (|188|)  that  as  n  -A  oo,  the  gravitational  coupling  g(n),  to 
this  order  in  the  Bessel  expansion,  has  to  scale  like 


a{n) 


1 

7^  ’ 


(191) 


so  that  the  expression  for  <  Atot  >  scales  like  n  or  NA,  with  the  expression  for  <  Aa  >  staying 
finite. 

The  result  of  Eq.  (119011  for  xa  then  implies 

1 


XA 


n 


9  Vn 


(192) 


in  the  same  limit  n  -A  oo.  In  view  of  Eqs.  (1 1871)  and  (11851)  with  n  ~  NA  ~  L2,  this  would  imply 
2/v  —  3  =  0,  and  thus  for  the  universal  critical  exponent  v  itself  v  =  |  =  0.666  to  first  order 
(m  =  1)  in  the  Bessel  function  expansion  of  Eq.  (116111  and  v  =  =  0.588  to  the  next  order 

(m  =  2)  in  the  same  expansion. 

With  some  additional  work  one  can,  in  fact,  completely  determine  the  asymptotic  behavior  of 
various  averages  for  large  /3  (small  g)  and  large  n.  First  one  notes  that  when  m  Bessel  functions 
are  included  in  the  expansion  for  the  wave  function  given  in  Eq.  (|161l).  beyond  the  leading  order 
one  at  strong  coupling,  one  obtains  a  wave  function  which  contains  powers  of  /3  up  to  /3m.  For  a 
given  fixed  m  one  then  finds  for  the  average  area  per  triangle  the  following  asymptotic  result 

1 


<  Aa> 


m+1  5 

gdm  1  2 


(193) 


up  to  terms  which  contain  higher  powers  of  1/n  (making  these  less  relevant  in  the  limit  n  -A  oo), 
and  also  up  to  terms  which  are  less  singular  in  g  for  small  g.  The  requirement  that  the  average 
area  per  triangle  be  finite  as  n  -A  oo  then  requires  that  the  coupling  g  itself  should  scale  with  n 
according  to 

S{n)  ~  — ^+r~  ■  (194) 


m+1 

TL  2(3m-l) 


For  the  area  fluctuation  itself  one  then  computes  in  the  same  limit 
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XA 
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3m- 2  n ^ 


(195) 
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again  to  leading  order  in  1  jn  and  1/g.  The  requirement  that  g(n )  scale  according  to  Eq.  (I194|)  then 
implies  from  Eq.  (11951)  that  the  area  fluctuation  diverges  in  the  limit  n  — >■  oo  as 


Xa(ji)  ~  n3m-! 


(196) 


By  comparing  with  Eqs.  (I184D  and  (11851)  one  obtains  immediately  for  the  exponent 

2  m  —  2 


a 

v 


3m  —  1 

and,  therefore,  from  the  scaling  relation  a.  =  2  —  3  v  finally 

6m  —  2 

v  = 


(197) 


(198) 


11  m  —  5 

One  can  now  take  the  limit  m  — >•  oo  [infinite  number  of  Bessel  functions  retained  in  the  expansion  of 
Eq.  (11611)] ,  which  leads  to  the  exact  result  for  the  correlation  length  exponent  of  2  +  1  dimensional 
quantum  gravity 


v  =  —  =  0.5454. 
11 


(199) 


The  derivation  shows  that  the  exponent  v  does  not  seem  to  depend  on  the  Euler  characteristic  x 


and,  therefore,  on  the  boundary  conditions 



i  Furthermore  one  can  compare  the  above  value  for 
v  with  the  (numerically  exact)  Euclidean  three-dimensional  quantum  gravity  result  obtained  over 
twenty  years  ago  in  [21],  namely  ^  ~  0.59(2).  It  would,  of  course,  be  of  great  interest  to  repeat  the 
above  Euclidean  lattice  calculation  in  order  to  refine  the  estimate  and  improve  on  the  statistical  and 
systematic  uncertainty.  The  exponent  v  is  expected  to  represent  a  universal  quantity,  independent 
of  short-distance  regularization  details  and,  therefore,  characteristic  of  gravity’s  universal  scaling 
properties  on  distances  much  larger  than  the  lattice  cutoff.  As  such,  it  should  apply  equally  to  both 
the  Lorentzian  and  the  Euclidean  formulation,  and  our  results  are  consistent  with  this  conclusion. 
Moreover,  in  3  +  1  dimensions  the  exponent  v  is  a  key  physical  quantity  as  it  determines  the  power 
for  the  running  of  the  gravitational  constant  G  [31]  and  for  the  Euclidean  theory  it  is  known  m 
that  the  universal  scaling  exponent  is  consistent  with  v  =  1/3. 

It  is  perhaps  worthwhile  at  this  point  to  compare  with  other  attempts  at  determining  the  critical 
exponent  v  in  three-dimensional  gravity.  The  latest  and  best  results  for  quantum  gravity  in  the 
perturbative  diagrammatic  2  +  e  continuum  expansion  using  the  background  held  method 
give  in  d  =  3  (e  =  1  and  central  charge  c  =  1) 


v-1  =  1  +  7  + 
5 


(200) 


8One  might  wonder  if  the  value  for  v  is  affected  by  the  choice  of  normalization  in  Eqs.  m  and  (11551) .  It  is  easy  to 
check  that  at  least  the  inclusion  of  a  weight  factor  Am ,  with  m  integer,  does  not  change  the  result  given  in  Eq.  Ifl99)l. 
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to  two-loop  order  and,  therefore,  v  ~  0.625,  with  a  substantial  uncertainty  of  about  fifty  percent 
(which  can  be  estimated  for  example  by  comparing  the  one-  and  two-loop  results).  On  the  other 
hand,  truncated  renormalization  group  calculations  for  gravity  directly  in  three  dimensions  [37,  [38] 
give  to  lowest  order  in  the  truncation  (i.e.  with  the  inclusion  of  the  cosmological  and  Einstein- 
Hilbert  terms  only)  the  estimate 


_  2d{d-2) 
~  d  +  2 


(201) 


and,  therefore,  in  d  =  3  the  value  v  ~  0.833.  This  last  result  is  also  affected  by  a  rather  substantial 
uncertainty  (again  as  much  as  fifty  percent),  which  can  be  estimated,  for  example,  by  including 
curvature-squared  terms  in  the  truncated  expansion.  Nevertheless,  and  in  light  of  the  uncertainties 
associated  with  the  various  methods,  it  is  very  encouraging  to  note  that  widely  different  calculations 
(on  the  lattice  and  in  the  continuum)  give  values  for  the  universal  scaling  exponent  v  that  are 
roughly  in  the  same  ballpark. 

From  Eq.  (11991)  one  obtains  the  fractal  dimension  for  a  gravitational  path  in  2  +  1  dimensions 


v  1  =  dp 


11 

IF 


1.8333... 


(202) 


This  is  slightly  smaller  than  the  value  for  a  free  scalar  field  dp  =  2,  corresponding  to  the  Brownian 
motion  (or  Wiener  path)  value.  It  is  closer  to  the  value  expected  for  a  dilute  branched  polymer  in 
the  same  dimension  [391 140],  and  the  best  match  at  this  point  seems  to  be  the  0(n)  vector  model 
for  n  =  —  1.  The  exact  value  v  =  6/11  for  2  +  1  gravity  would  then  suggest  a  connection  between 
the  ground  state  properties  of  quantum  gravity  and  the  geometry  of  dilute  branched  polymers  in 
the  same  dimension. 

In  light  of  the  results  obtained  so  far  it  is  possible  to  make  a  number  of  additional  observations. 
First,  note  from  Eq.  (11941)  that  as  n  — >•  oo,  the  critical  point  (or  renormalization  group  ultraviolet 
fixed  point)  moves  to  g  =  0 

9(n)  ~  -J7e  ■  (203) 

For  comparison,  a  variational  calculation  based  on  correlated  product  (Slater-Jastrow)  wave  func¬ 
tions  [1]  in  2  +  1  dimensions  gave 


9c  = 


4tt  x 


(204) 


N A  xAoOd  -  2)  ’ 

where  <to  >  2  was  a  parameter  associated  there  with  the  choice  of  functional  measure  over  edges. 
The  variational  result  of  Eq.  (12041)  can  be  compared  directly  with  the  result  of  Eqs.  (I165D  and  (I203D. 
for  x  =  2  and  Na  =  2n  +  2.  Thus  in  both  treatments  the  limiting  value  for  the  critical  point  for  g 
in  2  +  1  dimensions  is  zero,  gc  — >•  0  as  the  number  of  triangles  Na 


00. 
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Physically,  this  last  result  implies  that  there  is  no  weak  coupling  phase  ( g  <  gc,  or  in  terms 
of  Newton’s  constant  G  <  Gc ):  the  only  surviving  phase  for  gravity  in  three  dimensions  is  the 
strongly  coupled  one  ( g  >  gc  or  G  >  Gc ).  Furthermore,  the  correlation  length  £  of  Eq.  (I183D  is 
finite  for  g  >  0  and  diverges  at  g  =  0.  In  particular,  the  weak  field  expansion,  which  assumes  g 
small,  is  expected  to  have  zero  radius  of  convergence.  In  a  sense  this  is  a  welcome  result,  as  in 
the  Euclidean  theory  the  weak  coupling  phase  was  found  to  be  pathological  and  thus  physically 
unacceptable  in  both  three  [23]  and  four  dimensions  mm-  It  would  seem,  therefore,  that  the 
Euclidean  and  Lorentzian  lattice  results  are  ultimately  completely  consistent:  quantum  gravity 
in  2  +  1  dimensions  always  resides  in  the  strong  coupling,  gravitational  antiscreening  phase;  the 
weak  coupling,  gravitational  screening  phase  is  physically  excluded.  In  addition,  the  exact  value 
for  v  determines,  through  standard  renormalization  group  arguments,  the  scale  dependence  of  the 
gravitational  coupling  in  the  vicinity  of  the  ultraviolet  fixed  point  [34]. 


10  Summary  and  Conclusions 

In  this  paper  we  have  discussed  the  form  of  the  gravitational  wave  function  that  arises  as  a  solution 
of  the  lattice  Wheeler-DeWitt  equation  [Eqs.  (1201) .  (12TT)  and  ([541)]  for  finite  lattices.  The  main  result 
was  the  wave  function  T  given  in  Eqs.  (11351).  (11471)  and  (I153D  with  strong  coupling  limit  (curvature 
term  absent)  corresponding  to  the  choice  of  parameter  (5=0. 

To  summarize,  and  for  the  purpose  of  the  following  discussion,  the  wave  function  T  given  in 
Eq.  ([153D  can  be  written  in  the  most  general  form  as 

T  ~  e~ix  iFi  (a,  b,  2ix)  (205) 

up  to  an  overall  normalization  constant  N,  and  with  parameters  related  to  various  geometric 

9  These  circumstances  are  perhaps  unfamiliar  in  the  gravity  context,  but  are  nevertheless  rather  similar  to  what 
happens  in  gauge  theories,  including  compact  Quantum  Electrodynamics  in  2  +  1  dimensions  m-  There,  the  theory 
always  resides  in  the  strong  coupling  or  disordered  phase,  with  a  finite  correlation  length  which  eventually  diverges 
at  zero  charge. 

10  Specifically,  the  universal  exponent  v  is  related  to  the  behavior  of  the  Callan-Symanzik  beta  function  for  Newton’s 
constant  G  in  the  vicinity  of  the  ultraviolet  fixed  point  by  (}'  (G)\g=gc  =  —  1/v.  Integration  of  the  renormalization 
group  equations  for  G  then  determines  the  scale  dependence  of  G(/x)  or  G(D)  in  the  vicinity  of  the  ultraviolet  fixed 
point.  Concretely,  v  determines  the  exponent  in  the  running  of  G.  One  finds  G(D)  ~  (£2HH)  1^2i/,  with  □  =  grM"VMV„ 
the  covariant  d’Alembertian  and  £  the  renormalization  group  invariant  correlation  length.  A  broader  discussion  of 
renormalization  group  methods  as  they  apply  to  quantum  gravity  can  be  found  for  example  in  [19] . 
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invariants 


a  =  \  Na  —  — 'jzz — -  \  =  4  Na  —  — j= —  [  d2y  Ng  R 

4  V\G  4  2V2XG  J  V 

b  =  \N A 

x  =  ^=rAtot  =  J  d2y  y/g  .  (206) 


In  the  above  definitions  one  can  trade,  if  one  so  desires,  the  total  number  of  triangles  Na  for  the 
total  area 


Na  = 


1 


<Aa>  Atot  ~  <Aa 



(207) 


Use  has  been  made  of  the  relationship  between  various  coupling  constants  (<?,  G,  (3,  A,  A)  to  reexpress 
the  wave  function  ^  in  slightly  more  general  terms,  as  a  function  of  the  original  couplings  A  and  G 
appearing  in  the  original  form  of  the  Wheeler-DeWitt  equation  [see  for  example  Eqs.  (1421)  and  (1441)] . 
We  did  show  that  an  equivalent  form  for  the  wave  function  'I'  can  be  given  in  terms  of  Coulomb 
wave  functions  [see  Eq.  (11471)].  with  argument 

\/2  7T  X  1 


p  = 


Vx  G  2V2XG 


d~x  y/g  R 


(208) 


and  x  defined  as  in  Eq.  (I206D. 

The  above  wave  function  is  exact  in  the  limit  of  large  areas  and  completely  independent  of  the 
weak  field  expansion.  Nevertheless  it  is  only  correct  to  some  low  order  in  the  same  expansion  in  the 
limit  of  small  areas.  This  situation  was  interpreted  as  follows.  For  large  areas  one  has  a  very  large 
number  of  triangles,  and  the  short  distance  details  of  the  lattice  setup  play  a  vanishingly  small  role 
in  this  limit.  One  recognizes  this  limit  as  being  relevant  for  universal  scaling  properties,  including 
critical  exponents.  For  small  areas  on  the  other  hand  a  certain  sensitivity  to  the  short  distance 
properties  of  the  lattice  regularization  persists,  and  thus  a  universal  behavior  is,  not  unexpectedly, 
hard  to  achieve.  In  any  case  this  last  limit,  in  the  absence  of  a  truly  fundamental  and  explicit 
microscopic  theory,  is  always  expected  to  be  affected  by  short  distance  details  of  the  regularization, 
no  matter  what  its  ultimate  nature  might  be  (a  lattice  of  some  sort,  dimensional  regularization,  or 
an  invariant  continuum  momentum  cutoff,  etc.) 

In  principle  any  well-defined  diffeomorphism-invariant  average  can  be  computed  using  the  above 
wave  functions.  This  will  involve  at  some  point  the  evaluation  of  a  vacuum  expectation  value  of 
some  operator  O(g) 

I dg[g]  O(gij)  \^[gij]\2 
/  dg[g\ 


(*\6(gm 
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where  dfi[g\  is  the  appropriate  functional  measure  over  three- metrics  g^j.  Evaluating  such  an  average 
is,  in  general,  non-trivial,  as  it  requires  the  computation  of  a  (Euclidean)  lattice  path  integral  in 
one  dimension  less 

<tf|6(5)|tf)  =  A r  J  dg[g\d(gij )  exp  {-Seff[g]}  (210) 

with  Seff[g\  =  —  In  |'I'[(7jj]|2  and  A f  a  normalization  constant.  The  operator  O(g)  itself  can  be 
local,  or  nonlocal  as  in  the  case  of  the  gravitational  Wilson  loop  discussed  in  |42| .  Note  that  the 
statistical  weights  have  many  zeros  corresponding  to  the  nodes  of  the  wave  function  T,  and  that 
Seff  is  infinite  there. 

In  the  previous  sections  we  have  shown  that  the  wave  function  allows  one  to  calculate  a  number 
of  useful  and  physically  relevant  averages  and  fluctuations,  which  were  later  extrapolated  to  the 
infinite  volume  limit  of  infinitely  many  triangles.  It  was  found  that  these  diffeomorphism-invariant 
observables  point  in  2  +  1  dimensions  to  the  existence  of  a  fixed  point  (a  phase  transition  in  statis¬ 
tical  held  theory  language)  at  the  origin,  Gc  =  0.  One  concludes,  therefore,  that  the  weak  coupling 
(or  gravitational  screening)  phase  has  completely  disappeared  in  the  lattice  nonperturbative  formu¬ 
lation  and  that  the  theory  resides  in  the  strong  coupling  phase  only.  By  contrast,  in  the  Euclidean 
theory  it  was  found  in  [23]  that  the  weak  coupling  or  gravitational  screening  phase  exists  but  is 
pathological,  corresponding  to  a  degenerate  branched  polymer.  A  similar  set  of  results  is  found  in 
the  four-dimensional  Euclidean  theory,  where  the  weak  coupling,  gravitational  screening  phase  also 
describes  a  branched  polymer.  []_] 

The  calculations  presented  in  this  paper  and  in  |Tj  can  be  regarded,  therefore,  as  consistent 
with  the  conclusions  reached  earlier  from  the  Euclidean  framework,  and  no  new  surprises  arise 
when  considering  the  Lorentzian  2  +  1  theory.  Furthermore,  we  have  emphasized  before  that  the 
results  obtained  point  at  a  nonanalyticity  in  the  coupling  at  G  =  0,  signaling  a  strong  vacuum 
instability  of  quantum  gravitation  in  this  dimension.  In  view  of  these  results  it  is  therefore  not 
surprising  that  calculations  that  rely  on  the  weak  field,  semiclassical  or  small  G  expansion  run  into 
serious  trouble  and  uncontrollable  divergences  very  early  on.  Such  an  expansion  does  not  seem  to 
exist  if  the  non-perturbative  lattice  results  presented  here  are  taken  seriously.  The  correct  physical 
vacuum  apparently  cannot  in  any  way  be  obtained  as  a  small  perturbation  of  flat  or  near-flat 
spacetime. 

Let  us  add  here  a  few  final  comments  aimed  at  placing  the  present  work  in  the  context  of 
previous  calculations  for  the  same  theory.  A  number  of  attempts  have  been  made  over  the  years 


11  The  nature  of  solutions  to  the  lattice  Wheeler- DeWitt  equation  in  3  + 1  dimensions  will  be  discussed  in  a  separate 
publication  m- 
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to  obtain  an  estimate  for  the  gravitational  wave  functional  \h[g]  in  the  absence  of  sources.  These 
generally  have  relied  on  the  weak  field  expansion  in  the  continuum,  as  originally  done  in  [8j  [9]. 
Thus,  for  example,  one  finds  in  3  +  1  dimensions 


y[hTT}  =  Af  exp  | —  l  J  d3k  k  hJ^F (k)  ,  (211) 

where  tif^ (k)  is  the  Fourier  amplitude  of  transverse-traceless  modes  for  the  linearized  gravitational 
field  in  four  dimensions.  The  above  wave  functional  describes  a  collection  of  harmonic  oscillator 
wave  functions,  one  for  each  of  the  infinitely  many  physical  modes  of  the  linearized  gravitational 
field.  As  in  the  case  of  the  electromagnetic  field,  the  ground  state  wave  functional  can  be  expressed 
equivalently  in  terms  of  first  derivatives  of  the  field  potentials  (the  corresponding  B  field  for  gravity), 
without  having  to  mention  Fourier  amplitudes,  as 


y[hTT] 
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d3x  /  d3y 
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(212) 


Clearly  both  of  the  above  expressions  represent  only  the  leading  term  in  an  expansion  involving 
infinitely  many  terms  in  the  metric  fluctuation  hij  (and  since  they  apply  to  an  expansion  about  flat 
space,  the  cosmological  constant  term  does  not  appear  either).  Now,  in  2  +  1  dimensions  the  above 
expressions  become  meaningless,  since  there  cannot  be  any  transverse-traceless  modes.  The  only 
expectation  that  remains  true  is  that  the  wave  functional  should  still  depend  on  physical  degrees 
of  freedom  only:  it  should  be  a  function  of  the  intrinsic  geometry  of  3-space  and  should  not  change 
under  a  general  coordinate  change. 

If  one  restricts  oneself  to  local  terms  a  number  of  invariants  are  possible  in  2  +  1  dimensions. 
In  principle,  the  wave  function  could  depend  on,  besides  the  total  area 


Atot 
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and  curvature 
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other  invariants  such  as 

rn  =  J  d2x  yfg  Rn 


(215) 


with  n  an  integer.  The  latter  result  follows  from  the  fact  that  in  d  =  2,  both  the  Riemann  and 
Ricci  tensors  only  have  one  component,  related  to  the  scalar  curvature, 
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Nonlocal  terms  are  possible  as  well  involving  inverse  powers  of  the  covariant  d’Alembertian  □,  but 
these  do  not  seem  to  play  a  significant  role  in  the  lattice  theory. 

Now,  the  relevant  Euclidean  theory  for  the  present  work  is,  of  course,  gravity  in  three  (2  +  1) 
dimensions.  But  in  three  dimensions  the  Riemann  and  Ricci  tensor  have  the  same  number  of 
algebraically  independent  components  (6)  and  are  related  to  each  other  by 

=  e'"'  (R\  -  R'.)  (217) 

The  field  equations  then  imply,  using  Eq.  (12171) ,  that  the  Riemann  tensor  is  completely  determined 
by  the  matter  distribution  implicit  in  T^y, 

Rpupu  —  8 TtG  fjjip  Tl/(J  +  fji/a  T//p  T  Spa  Tup  Sup  Tp.rr  T  T  {Spcr  Sup  Spf>  Suer)]  (218) 


In  empty  space  =  0,  which  then  implies  for  zero  cosmological  constant  the  vanishing  of  Riemann 
there 

Rpupa  =  0  .  (219) 


As  a  result  in  three  dimensions  classical  spacetime  is  locally  flat  everywhere  outside  a  source,  gravi¬ 
tational  fields  do  not  propagate  outside  matter,  and  two  bodies  cannot  experience  any  gravitational 
force:  they  move  uniformly  on  straight  lines.  There  cannot  be  any  gravitational  waves  either:  the 
Weyl  tensor,  which  carries  information  about  gravitational  fields  not  determined  locally  by  matter, 
vanishes  identically  in  three  dimensions. 

What  seems  rather  puzzling  at  first  is  that  the  Newtonian  theory  seems  to  make  perfect  sense 
in  d  =  3.  This  can  only  mean  that  the  Newtonian  theory  is  not  recovered  in  the  weak  field  limit  of 
the  relativistic  theory.  To  see  this  explicitly,  it  is  sufficient  to  consider  the  trace-reversed  form  of 
the  field  equations, 

Rpu  =  8? tG  ^  T/w  -  T  \ 

with  T  =  T\.  in  the  weak  field  limit.  In  the  linearized  theory,  with  h^u  =  Sp,u  ~  Vpu ,  and  in  the 
gauge  Va hx^  —  ^Vflhxx  =  0,  one  obtains  the  wave  equation 
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(221) 


with  T/,y  the  linearized  stress  tensor.  After  neglecting  the  spatial  components  of  rfW  in  comparison 
to  the  mass  density  too,  and  assuming  that  the  fields  are  quasistatic,  one  obtains  a  Poisson  equation 
for  h00, 

d  3 

V2/ioo  =  —167 TG  - — -  tqo  (222) 


53 


In  four  dimensions  this  is  equivalent  to  Poisson’s  equation  for  the  Newtonian  theory  when  one 
identifies  the  metric  with  the  Newtonian  field  4>  in  the  usual  way  via  /ioo  =  —2 <f>.  But  in  three 
dimensions  such  a  correspondence  is  obstructed  by  the  fact  that,  from  Eq.  (12221).  the  nonrelativistic 
Newtonian  coupling  appearing  in  Poisson’s  equation  is  given  by 

^Newton  =  G  (223) 

and  the  mass  density  too  completely  decouples  from  the  gravitational  field  /loo-  As  a  result,  the 
linearized  theory  in  three  dimensions  fails  to  reproduce  the  Newtonian  theory. 

In  a  complementary  way  one  can  show  that  gravitational  waves  are  not  possible  either  in  the 
linearized  theory  in  three  dimensions.  Indeed  the  counting  of  physical  degrees  of  freedom  for  the 
d-dimensional  theory  goes  as  follows.  The  metric  g^iV  has  ^d(d  +  1)  independent  components;  the 
Bianchi  identity  and  the  coordinate  conditions  reduce  this  number  to  \d(d+  1)  —  d  —  d  =  \d(d—  3), 
which  gives  indeed  the  correct  number  of  physical  degrees  of  freedom  (two)  corresponding  to  a 
massless  spin  two  particle  in  d  =  4,  and  no  physical  degrees  of  freedom  in  d  =  3  (and  minus  one 
degree  of  freedom  in  d  =  2,  which  is  in  fact  incorrect).  Nevertheless,  investigations  of  quantum 
two-dimensional  gravity  have  uncovered  the  fact  that  there  can  be  surviving  degrees  of  freedom  in 
the  quantum  theory,  at  least  in  two  dimensions.  The  usual  treatment  of  two-dimensional  gravity 
[MJ  starts  from  the  metric  in  the  conformal  gauge  gtiv(x)  =  e^(x)g^ll/,  where  g ^  is  a  reference 
metric,  usually  taken  to  be  the  flat  one.  The  conformal  gauge-fixing  then  implies  a  nontrivial 
Faddeev-Popov  determinant,  which,  when  exponentiated,  results  in  an  effective  Liouville  action, 
with  a  potential  term  coming  from  the  cosmological  constant  contribution.  One  would,  therefore, 
conclude  from  this  example  that  gravity  in  the  functional  integral  representation  needs  a  careful 
treatment  of  the  conformal  degree  of  freedom,  since  in  general  its  dynamics  cannot  be  assumed 
to  be  trivial.  The  calculations  presented  in  this  paper  show  that  this  is  indeed  the  case  in  2  +  1 
dimensions  as  well. 
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